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CHAPTER  I 
INTRODUCTION 


Among  the  various  problems  of  electromagnetic  scattering  by  arbi¬ 
trary  bodies,  only  a  few  simple  shapes  have  led  to  exact  solutions. 

The  diffraction  of  an  electromagnetic  plane  wave  by  a  conducting  circu¬ 
lar  disk  long  remained  unsolved.  The  presence  of  an  edge  led  to  much 
greater  complications  than,  for  example,  the  sphere  which  was  solved  by 
Mie[l]  very  early  in  this  century. 

First,  some  limiting  cases  were  derived  for  the  disk  problem. 
Approximate  solutions  were  obtained  by  Kirchhoff  diffraction  theory 
for  disks  of  large  electrical  circumference  kd  where  k  is  the  wave 
number  and  d  the  diameter  of  the  disk  and  by  Lord  Rayleigh[2]  for 
disks  of  small  electrical  circumference.  After  the  second  world  war, 
extensive  work  was  directed  at  the  disk  problem,  and  two  complete  gen¬ 
eral  solutions  were  developed.  Meixner  and  Andrejewski [3]  used  Hertzian 
vector  potentials  to  compute  the  scattered  field.  They  derived  the  form 
of  the  potential  from  the  appropriate  boundary  conditions,  and  the  uni¬ 
queness  of  the  solution  was  insured  by  the  edge  condition  developed  by 
Meixner[4].  He  showed  that  the  energy  density  of  the  total  field  in  the 
vicinity  of  an  edge  must  be  integrable.  This  insists  that  the  compon¬ 
ents  of  the  fields  vary  at  most  as  s"*'2,  where  s  is  the  distance  to  the 
observation  point  from  the  edge.  Bouwkamp[5]  showed  that  the  tangential 
electric  field  is  zero  on  the  edge,  that  it  approaches  this  value  as 
s1'2,  and  that  the  tangential  magnetic  field  remains  finite.  Tangen¬ 
tial,  in  this  context,  means  parallel  to  the  edge.  This  solution  has 
been  checked  by  numerous  computations  made  in  recent  years.  Since  the 
development  by  Hodge[6]  of  an  efficient  method  of  computing  the  spher¬ 
oidal  eigenvalues,  even  more  efficient  computers  have  allowed  extensive 
comparison  between  experimental  data  and  computational  results  with  very 
good  agreement. 

Shortly  after  Meixner,  Flammer[7]  derived  another  solution  using 
oblate  spheroidal  vector  wave  functions  to  expand  the  fields.  As  in  the 
Mie  solution  for  the  sphere,  the  fields  of  the  incoming  plane  wave  were 
expanded  in  terms  of  the  vector  wave  functions.  The  boundary  conditions 
on  the  surface  of  the  disk  determined  the  scattering  components.  In 
order  to  insure  that  the  tangential  component  of  the  electric  field  is 
zero  on  the  edge,  two  sets  of  vector  wave  functions  were  used  and  their 
relative  weights  in  the  expansion  of  the  incident  plane  wave  were  deter¬ 
mined  by  the  edge  boundary  condition.  This  led  to  a  unique  solution. 
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The  purpose  of  this  work  is  to  compare  Meixner's  and  Flammer's 
solutions.  Apparently,  no  numerical  results  have  ever  been  obtained 
from  Flammer's  formal  solution.  Thus,  it  is  necessary  to  establish  that 
the  formal  solution  is,  in  fact,  valid.  If  this  validity  is  established, 
the  vector  wave  function  formalism  may  be  more  convenient  for  examina¬ 
tion  of  such  characteristics  as  the  surface  current  distributions  and 
the  natural  resonances.  The  reference  work  used  for  Meixner's  solution 
is  Hodge's  version.  Hodge[8]  used  the  notation  and  normalization  of  the 
spheroidal  functions  introduced  by  Flammer.  In  order  to  match  Hodge's 
geometry,  Flammer's  solution  is  rederived  for  the  same  configuration. 

The  bistatic  normal  incidence  case  is  more  deeply  studied,  and  problems 
encountered  in  the  derivation  of  the  solution  are  pointed  out.  Compu¬ 
tations  made  in  the  normal  incidence  case,  as  studied  by  Flammer,  are 
compared  to  the  results  obtained  by  Hodge[9].  In  addition,  Meixner's 
fields  are  expressed  in  terms  of  spheroidal  vector  wave  functions  as 
well  as  the  vector  potentials.  Another  proof  of  Meixner's  solution  with 
vector  wave  functions  is  given. 
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CHAPTER  II 

SPECIFICATIONS  OF  THE  PROBLEM 


A.  Geometry 

The  scattering  object  is  an  infinitely  thin,  perfectly  conducting 
circular  disk  of  radius  a.  It  lies  in  the  x-y  plane  of  a  right-handed 
Cartesian  coordinate  system  and  is  centered  at  the  origin. 

The  following  study  will  be  made  in  the  oblate  spheroidal  coordi¬ 
nate  system  (n,C>4>),  Figure  2.1,  in  which  the  disk  is  the  surface  repre¬ 
sented  by  c-0. 


*7  =  0 


1  <0 


Figure  2.1.  The  oblate  coordinate  system. 


The  transformation  from  oblate  spheroidal  to  Cartesian  coordinates 
are,  as  given  by  Flammer[10] , 
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X 


'j-2a 


=  a[(l-n2)0+f2)]1/2  cos* 
jv  -  a[(l-n2)(lH2)]1/2  s  i  n^>  (1) 

[2  =  a  nC 


In  Appendix  A,  the  transformations  between  the  different  coordi¬ 
nate  systems  -  Cartesian,  cylindrical  and  spherical  -  and  the  oblate 
spheroidal  coordinate  system  are  summarized.  Their  limiting  forms  are 
listed  for  three  cases:  at  large  distances  from  the  disk  -£  large  -  on 
the  surface  of  the  disk  and  in  the  vicinity  of  the  edge. 


B.  Problem 


Because  of  the  rotational  symmetry  about  the  z-axis  it  is  possible, 
without  loss  of  generality,  to  choose  an  arbitrary  plane  of  incidence. 

We  will  use  the  notation  of  Hodge's  work[9].  The  incident  plane  wave 
direction  lies  in  the  x-z  plane,  coming  from  the  x>0  half  plane,  Figure 
2.2.  flammer  used  the  y-z  plane  as  plane  of  incidence,  instead.  The 
whole  system  is  located  in  free-space. 


Figure  2.2.  Geometry. 
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We  will  use  the  following  symbols: 
c=ka  -  electrical  circumference 
e0  -  angle  of  incidence 
a  -  polarization  angle 
e1u)t  -  time  dependence 
The  incident  wave  vector  F1  is  expressed  by: 


F1  =  k[-sineQ  ex  -  cose0  ez]  (2) 


The  polarization  of  the  incident  electric  field  is  measured  bj 
the  angle  a  between  the  plane  of  incidence  and  the  incident  electric 
field,  E i ,  in  such  a  way  that  a  is  the  conventional  azimuthal  angle, 
of  a  coordinate  system  defined  by  the  unit  vectors  (cose0  ex-sine0  ez, 
iy,  sine0  ex+cose0  ez)  as  shown  in  Figure  2.3.  Parallel  polarization 
is  obtained  for  a=0  and  perpendicular  polarization  for  a=Tr/2.  The 
incident  electric  field  is  given  by 


Ep  =  Eo(cos0o  cosa  ex  +  sina  iy  -  sin90  cosa  ez) 

•  e'1^1  ‘ r  ela,1; 


(3) 


where  r  is  the  usual  position  vector. 

Note  that  Flammer[7]  characterizes  the  incident  wave  by  the  angle 
r  between  the  positive  direction  of  the  propagation  vector  and  the  posi¬ 
tive  z  axis  instead  of  the  angle  of  incidence  90.  Those  angles  are  re¬ 
lated  as  follows:  c=ir-e0. 

In  the  spheroidal  coordinate  systems,  the  boundary  conditions  sat¬ 
isfied  by  the  electromagnetic  field  on  the  surface  of  the  disk  and  at 
the  edge  are  on  the  disk,  £=0: 

(r+P)  •  en  =  0 

(r+P)  •  e^  =  0  (4) 

(H1 +H3)  •  er  =  0 

where  E3  is  the  scattered  electromagnetic  field. 
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Figure  2.3.  Polarization 
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At  the  edge,  c= 0  and  a=0,  the  energy  density  of  the  total  field 
must  be  integrable.  The  components  of  the  fields  may  become  infinite 
at  most  as  s'1/2,  where  s  is  the  distance  from  the  edge  to  the  observa¬ 
tion  point.  Bouwhamp[10]  showed  the  following  more  restrictive  condi¬ 
tions  on  the  ^-component  of  the  fields: 


(r+£s)  •  e^  =  0(s1/2) 

OT+H5)  •  e+  =  0(1)  (5) 

Those  orders  of  variation  as  function  of  s  are  valid  only  in  the 
vicinity  of  the  edge,  s<<a,  for  the  curved  edge  of  a  disk. 

In  this  work,  we  will  refer  to  these  behaviors  of  the  components 
of  the  total  electromagnetic  field  in  the  vicinity  of  the  edge  as  the 
edge  condition.  The  ^-components  will  also  be  refered  to  as  tangential 
to  the  edge,  or  parallel  to  the  rim. 

We  will  see  that,  in  Flammer's  solution,  the  condition  on  [E1+Er]({j 
is  equivalent  to  the  whole  edge  condition.  This  is  due  to  the  particular 
choice  of  the  vector  wave  functions  used.  This  will  not  be  true  for 
Meixner's  solution,  where  all  the  conditions  on  the  components  of  the 
field  will  be  needed. 

Before  we  present  Flammer's  and  Meixner's  solution  we  need  to 
introduce  the  functions  that  we  are  going  to  work  with  and  some  of  their 
properties.  In  Chapter  III  we  will  define  the  spheroidal  functions,  and 
in  Chapter  IV  we  will  derive  various  plane  wave  expansions.  In  the  fol¬ 
lowing  work,  the  time  dependence  eiwt  will  be  omitted. 
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CHAPTER  III 

THE  SPHEROIDAL  FUNCTIONS 


In  this  chapter  the  scalar  wave  function  solution  of  the  scalar 
wave  equation  will  be  introduced  and  some  properties  of  the  spheroidal 
radial  and  angular  functions  will  be  summarized.  The  large  radial  argu¬ 
ment,  c>,,  approximation  of  the  vector  wave  function  solutions  of  the 
vector  wave  equation  will  be  listed. 


A .  Differential  Equation  and 

Scalar  Wave  Function 


In  the  oblate  spheroidal  coordinate  system,  the  scalar  wave  equa¬ 
tion  can  be  solved  by  the  method  of  separation  of  variables.  In  Appendix 
B  the  different  differential  operators  have  been  expanded  in  spheroidal 
coordinates.  The  scalar  wave  equation  is: 


V2.p+k  2<p  =  0 

1 


2n  +  l-n2  ijL  +  JS _ £1  +  1+C2 

f2+n2  3n  C2+n2  8n2  C2+n2  n2+C?  9£2 


4l- 


k2ip  =  0 


(1+L  )(1-n  )  3^j 

The  eigenfunctions  associated  with  the  eigenvalue  X|T)n  are: 


(6) 


</'  ( i )  ( n ,  C ,  <P ) 

mn 

o 


mn 


(~ic,n) 


r(  1 ) 

mn 


Hc,U) 


where  Smn(-ic,n)  is  the  angular  function 


S-nHc.r,)  -  V.  d;n(-ic)  p:+r(n) 


mn 


r=0,l 


m+rv 


as  defined  by  Flammer[10]. 


(7) 


(8) 
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Pm  (n)  is  the  Legendre  associated  function  as  used  in  spherical 
coordinates.  The  expansion  coefficients  dlJ!n(-ic)  are  computed  from  a 
recursion  relationship  as  seen  in  Flammer's  book[10]. 


The  prime  over  the  summation  is  a  notation  that  will  be  used 
throughout  the  following  work  to  denote  that  the  summation  must  involve 
only  the  even  or  odd  indexes  when  n-m  is  even  or  odd,  respectively. 

-  R^^(-ic.ic)  is  the  radial  function 

From  Flammer's  book[10]> 


(-ic.ic) 


(]n2)m 

£m  [  dmn(-ic) 
r=0,l  r 


(2m+r) ! 
r! 


co  1 

l  ir+ftw1  d"!n(-ic) 
r=0,l 

z(1)(cp  (9) 
r!  n+r 


The  index  (i)  denotes  the  kind  of  spherical  Bessel  function,  z^1). 
used:  i  =  l  for  jn,  i=2  for  nn,  i=3  for  hj,'  >  and  i=4  for  h^).  We  will 
refer  to  i  as  the  index  of  the  radial  function  in  the  rest  of  this  work. 
The  normalization  of  the  radial  and  angular  spheroidal  functions  are 
those  used  by  Flammer. 

In  the  following  two  sections,  we  will  summarize  some  properties 
of  the  angular  and  radial  spheroidal  functions. 


B.  The  Angular  Function  Smn(-ic,n) 


The  equation  satisfied  by  the  angular  functions  is: 


d_ 

dn 


[<’-',2)  Smn<-ic-n>]  *  (Amn«2"Z  ‘  ~t)  W-,'c-’>>  *  0 


-1  <  n  <  1 


The  angular  functions  are  orthogonal  with  weight  function  1: 

rl 


-1 


^mn ( - i c » n ) SmD ( - i c , n ) d  n  ^nD 


Jmp 


np 


(10) 


(ID 


where 
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(cTHc))2 


n  =  o  y  { r+2m) !  _ 1 _ 

mn  r=0  1  r!  2r+2m+l 


1 .  Evaluation  at  n=0 

From  Formula  (8),  or  more  easily  from  the  expansion  of  Smn(-ic,n) 
in  powers  of  ( 1  - n2 ) »  -  Formula  (3-2-7)  in  Reference  [10]  we  can  see 
that 

for  n-m  even  Smil(-ic,n)  is  an  even  function  of  n 

for  n-m  odd  Smn(-ic,n)  is  an  odd  function  of  n 


From  the  differential  equation  and  its  derivatives,  we  can  evalu¬ 
ate  the  derivatives  of  Smn(-ic,n)  for  zero  argument: 


T  CW'ic.O)]  =  -  (Amn-m2)  Smn(-ic,0) 


3  (Smn(-ic,0))  -  (Amn-m  -2)  Smn(-ic,0) 


7 r  )  '  t(Arnn-m^)(x[nn-m2-6)  +  2m  -2c  ]  Smn(-ic,0) 


Thus,  the  evaluation  of  Smn(-ic,n)  and  its  derivatives  around  n=0  can 
then  be  described  by  taking  into  account  the  odd  and  even  properties: 
n  SMALL  n«l 


2 

SmnHc,n)  =  Smn(-ic,0)  Cl  -  Umn-m2)]  +  0(n4) 


S;n(-ic,n)  =  Smn(-ic,0)  [-  n(*mn-m2)]  +  0(n3) 


Smn(_ic,n) 


-  Smn(“ic»0)  [-( Amn_l^2)+  9  ^^mn”*'1  ) ( A^^~m  -6) 


+  2m2-2c2]]  +  0(n4) 
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where  the  normalization  has  defined  Smn(-ic,0)  =  P™(0). 
n-m  odd 

/ 

Snm(-ic»n)  =  S^,n(-ic ,0)[n]  +  0  (n3) 

<s^n(-ic>n)  «  s;n(-ic,0)  n  -  £.  (xmn-m2-2)]  +  0(n4)  (14) 

smn(-ic.n)  =  s^n(-ic ,0)  [-n(Amn-m2-2)]  +  0(n3) 
where  the  normalization  has  defined  S^n(-ic,0)  =  ^[P™(0)] 

2.  Evaluation  at  n=±l 

The  power  series  expansion  is  led  by  (l-n2)m/^  due  to  the  P{{j+r(n) 
functions.  For  m  greater  than  zero. 


m  i  1  Smn(-ic,±l)  =  0 
m  >  3  S^n(-i c,±l )  =  0 


The  evaluation  for  m=0  and  1  has  been  given  by  FlammerpO]. 

S0„Hc,l)  *  V  «?nHO 

r=0,l 

S0n("ic’l )  =  \  1‘  ^(^+1)  d°n(-ic) 

L  r=0,l 

lim(l-n2)1/2  d  Csln(1c,n)3  *  -  i  V  (r+l)(r+2)  d^n(-ic) 
n+1  dr)  m  c  r=0,l  r 

Si  (-ic,n)  ”  7n 

lim  21/2  *  ^r  Pp+]0) 

n-*-l  (1-n  )  r=0,l 

1  l  (r+l)(r+2)din(-ic) 

2  r=0,l  r 


(15) 


(16) 


C.  The  Radial  Function  Rj(>)(-ic,iO 


The  equation  satisfied  by  the  radial  functions  is: 

2 


sr[('w-2)  at  ^’hc.ki]  -  1 


A  r2f2 
Amn  c  ^ 


in 


R^Hc.U)  =  0, 


o  <  e 


(17) 


Unlike  their  spherical  counterparts,  Bessel,  Neuman  and  Hankel 
functions  -  the  oblate  spheroidal  radial  functions  and  their  derivatives 
are  well  defined  over  the  whole  range  of  values  of  f;.  The  normalization 
coefficients  of  these  functions  have  been  defined  by  Flammer[10]  so  that 
they  match  convenient  formulas  for  large  values  of  the  argument  c£, 


(18) 


With  those  definitions,  we  see  that  the  suitable  functions  for 
the  series  expansion  of  an  incoming  plane  wave  will  involve  R^),  and 
outward  traveling  waves  will  involve  R^' . 


1 .  Wronskian 

The  Wronskian  of  the  differential  equation  is,  as  computed  by 
Flammer[10] : 


£><-10.10 


I^Hc.ic) 


dR^Hc.U) 

mn _ 

d£ 


=  + 


1 

c(lH2) 


(19) 


This  can  be  written  in  terms  of  the  R^,  and  R^  functions,  using 
the  property  ^^R^-iR^) ,  and  the  fact  that  the  Wronskian  of  two 
linearly  dependent  solutions  is  zero: 
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t 


tfHc.U)  ^  t^Mc.ic)]  -  ^Mc.U)  k  [R^Hc.ic)] 


c(l+C2) 


(20) 


2.  Evaluation  at  j=0 

In  his  book.  Flamer  lists  the  values  of  the  radial  functions  and 
their  derivatives  at  £=0.  Appendix  C  of  this  work  contains  some  correc¬ 
tions  to  those  formulas.  After  presenting  some  general  properties,  we 
focus  on  R^)  and  RW  separately. 

From  the  differential  equation  and  its  derivatives  we  can  evaluate 
the  derivatives  of  )  (-ic,  U)  at  £=0: 


z?  -  (v,-™2-2)  -He 


Sun 

dl&h-Ic.lO) 


dt, 

rR(i) 


tRmn,]("ic»i0)  =  C(Amn”m2)(Amn-m2-6)  +  2m2-2c2]  R^t-ic.iO) 

where  c=ka  (21) 


a.  Properties  of  R^ 

Flammer[10]  develops  a  power  series  expansion  of  R^P(-ic,i£). 
This  function,  if  extended  over  the  domain  -00<£<00>  would  lead  to  an 
even  or  odd  function  for  n-m  even  or  odd,  respectively.  We  can  also 
conclude  that,  for  k>0: 

~T  [Rj^H-ic.iO)  =  0  (22) 

for  k  even  or  odd  when  n-m  is  odd  or  even,  respectively. 
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In  the  vicinity  of  the  disk,  the  radial  functions  of  the  first  kind 
can  be  expanded  as  follows: 

r  Small  c^  <<  1 
n-m  even 

=  ^J)(-ic,iO)[l+(Amn-m2)  |ij  +  0(c4) 

4  Rmn),('ic’i0  =  ^^(-ic.iOJC^^n-m2)]  +  0(c3)  (23) 

Rii} "('ic,i°  =  Rmi)(-ic»i°)j Amn-m2  +  [(W"2'6)  +  2m2-2c2)] 

^  L  +  OU*)  J 


n-m  odd 


=  ^'(-ic.iOCc]  +  0(c3) 

Rmn} '  ("ic.ic)  =  R^} '  (-ic.iO)  ^l+(xmn-m2-2)|^  +  0(S4) 
(^^"(-ic.ic)  =  ^'(-ic.iO^Amn-m2^)  +  0(c3) 


b.  Properties  of  R^) 

From  Equation  (19),  the  Wronskian  of  rI^  and  R^4^  leads  us  to  the 
following  identities 


n-m  even  R^,4) '  (-ic,iO)  =  -  --/V  .  . 


c  C;(-ic,iO) 


n-m  odd  R^(-ic.iO)  =  - pryr1 - 

c  ii]n}  C-ic.iO) 


We  then  get  the  following  expansions  in  the  vicinity  of  the  disk: 
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For  c  Small ,  cf  <<  1 


n-m  even 


"  R^n^(-ic,iO)  |l  +  (\nn-rn2)j 


M7~.  r~  +  0(£3) 

c  ^mn  (-ic,iO) 


I  «tt>*C-ic.tO  -  -  — bi^w*2-*) 

+  ^W^^Hc.iO)  +  0(c3) 

^'‘Hc.U)  =  ^Hc.lO)  (26) 

[(Amn-m2)  +  ii  C(Amn-m2)(Amn-m2-6)  +  2m2-2c2)]J 

u  c  R^-ic.iO)  C(xmn-m2-2)  +  0(?3) 

Note  that  the  coefficients  of  R^-ic.iO)  in  these  equations  are 
(23).  t0  th0Se  of  (~ic,iO)  in  the  corresponding  equations  of  Formula 


n-m  odd 


R^(-ic,U)  =  — 7--|  \ 1 -  1  +  |- (Amn-m2)] 

c  (-ic.iO)  L  2  J 

+  5  Ri,S^(-ic,io)  +  OU3) 
Rmn)'  <-ic*i^)  *  ^  '  He, i0)  [  +  f  (W*2"2)] 


+  ~T0)‘!  •  4n,  S(Amn-m2)  +  0(C3) 

c  Km  He, 10) 


^  Hc,u>  ‘  navnco,  +  r 

♦  a»2-2c2]j  +  5  bAJ)  ,<-iC.iO){Amri-Inz-2}  *  0(i3) 
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Note  that  the  coefficients  of  (-ic,io)  in  these  equations 

are  equal  to  those  of  (-ic,iO)  in  the  corresponding  equations  of 
Formula  (24). 

3.  Evaluation  of  R^  for  large  arguments 
Equation  (18)  gives  the  large  argument  approximations 


knn^-ic.U) 


exp(-icf ) 


« [■•  •  a 


exp(-ic^) 


[Rir,nh("ic,ic)  - >  [c2  -  ~  -  ^  exPHcc) 

b 

In  the  far  field  expansions,  only  the  terms  with  dependence 
will  be  kept. 


D.  The  Vector  Wave  Functions 

Solenoidal  solutions  of  the  vector  wave  equation  can  be  expanded 
as  summations  of  spheroidal  vector  wave  functions.  Those  will  be  de¬ 
fined,  following  Flammer's  notation,  by: 


Mg'*' 1  '(n*C><f>) 


=  v  x  [ipg  ^ (n,C»cb)ea] 


Ng01^1  )(n,f  ,4>) 

omn 


I  v  x  M!(i)(n,C,<!>) 

K  £  mn 


where  e  is  either  a  constant  vector  in  Cartesian  coordinates  (a=x,y,z) 
or  the  position  vector  iNr  Sr(a=r). 

Extensive  work  has  been  done  on  these  functions  by  Flammer,  who 
lists  them  all  in  his  book[10].  In  contrast  with  the  orthogonality  of 
the  scalar  wave  functions,  the  spheroidal  vector  wave  functions  are  not 
orthogonal.  They  are  not  even  independent  as  will  be  seen  in  Chapter 
IV-C.  We  will  list  here  the  behavior  of  those  functions  on  the  surface 
of  the  disk  and  for  large  radial  arguments. 


1 .  Behavior  of  the  vector  wave 
functions  at  f,=Q 


We  have  seen  Chapter  III-C,  Equation  (21),  that,  on  the  surface 
of  the  disk,  all  radial  functions  and  their  derivatives  can  be  expressed 

in  terms  of  (-ic,iO)  and  R^i)  (-ic.iO).  We  can  show  that,  on  the 
surface  of  the  disk,  each  component  of  any  vector  wave  function  defined 
by  Equation  (29)  is  proportional  to  only  one  of  these  values.  We  list 
the  different  cases. 


a.  Tangential  components  at  g=0 


[M  x(i)]  ,  [fixM)]  ,  ,  [M^(i)]  ,[N"(i)]  and 

<ynn  n  emn  f  emn  i  emn  <P  emn  n 

v  0  0  0  0 

are 
r>mn 


/  i  \  » 

proportional  to  R^n;  (-ic,iO). 


(30) 


:;mn 


•  [SXe(i>]*  '  [Sf  )]n  •  •  [^’]n  • 

^mn  ^mn  *  ^mn  Qmn  v  ^  In  ' 

are  proportional  to  R^^ (-ic.iO) .  (31) 


b.  Normal  components  at  C=0 


[Nex(i)L 
wnn  s 

0 

,  [N/Oh 
> 

are  proportional  to  R^^  (-ic,iO) 

(32) 

[Nez^^]  are  proportional  to 

omn  ? 

„mn  ^  , 

i&’Mc.to) 

(33) 

We  notice  that  at  the  surface  of  the  disk  the  tangential  components 
of  any  of  these  vector  wave  functions,  and_the  normal  component  of  their 
curl  -  that  is,  for  example,  [N01]^  ^  and  [Ma]£  or  vice-versa  -  have  the 
same  radial  dependence.  That  ensures  that  both  boundary  conditions  on 
E  and  H  at  the  surface  of  the  disk  are  matched  by  the  same  expansion. 
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2.  Behavior  of  the  vector  wave  functions 
of  index  i=4  for  large  arguments  of 
the  radial  part 

In  this  section,  we  list  the  large  radial  argument  approximations 
of  the  vector  wave  functions,  keeping  only  the  1/f  terms.  From  Equation 
(29), 


fMeX(4)  =  f  f-  Smn(-ic>n){-sin4>  en+n  cos*  e^}  ™*m*  e‘1cC 
0mn 


Mey(4)  =  If”  Smn(-ic,n){cos<i»  en+n  sin*  e  l  ^m*  e-icz 
r,mn  ■ 


"eZ(4>  "  '  I  W--.n)  f  “J*  e-'« 

0mn 


(34) 


,-n-l 


cos*  en+sin*  e^}  ^m*  e"lC’ 


N|*‘4)  ■  jf-W-ic.nh, 
y(4)  *  2  i^J_  Smn(-ic,n){n  sin*  en-cos*  e^>  c^sm*  e~ic^ 


e  H  7 -  mn 

nmn  d  f 


sin 


z(4)  *  _  <•  /T- n2  Smn(-ic,n)  111  cosm*  e  e~icS 

mn  n  1,1  •  *  r  ~  A  n 


-mn 


S  sin 


We  now  have  enough  tools  to  find  the  expansions  of  scalar  and 
vector  plane  waves  in  terms  of  the  spheroidal  functions. 
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CHAPTER  IV 

PLANE  WAVE  EXPANSION 


In  order  to  expand  the  incident  fields,  we  need  the  expansions 
of  vector  plane  waves  in  terms  of  the  spheroidal  vector  wave  functions. 
First,  we  will  derive  the  scalar  plane  wave  expansion  in  terms  of  scalar 
wave  functions.  Taking  its  curl  associated  with  different  basis  vectors 
will  lead  us  to  various  vector  plane  wave  expansions.  The  special  cases 
of  normal  incidence  to  the  disk  are  reviewed. 


A.  Scalar  Plane  Wave 


By  scalar  plane  wave  we  denote  the  function  exp[ik(x  sinoQ  + 
z  cosGg)]  in  our  geometry.  We  do  not  know  enough  about  the  angular 
spheroidal  functions  to  directly  obtain  the  coefficients  of  the  expan¬ 
sion  of  the  scalar  plane  wave  in  terms  of  spheroidal  functions  by 
applying  orthogonality  on  the  trigonometric  and  angular  functions.  We 
will  use  another  method  as  derived  by  Flammer[10].  The  plane  wave  is 
considered  to  be  due  to  a  point  source  removed  to  infinity  in  the 
direction  of  arrival  of  the  incident  wave.  The  scalar  plane  wave  will 
then  be  obtained  from  the  asymptotic  form  of  the  free-space  Green's 
function. 

Because  of  our  choice  of  the  time  dependence,  e^ut,  the  Green's 
function  used  here  will  be  the  complex  conjugate  of  that  obtained  by 
Flammer[10] : 

exp(-jk|r' -r' I ) 

4ir|r'-r'  | 


(35) 


CO  03 


lk  l_  I  ITT1 


2tt 


m=0  n=m  ron 


Ri!i^(-ic,U'  )Rmn^(-ic,ic)  cosm( 4,-4. * )  when  £‘>C, 


where  ( n *  '  *4> ' )  are  the  coordinates  of  the  source  point. 

For  a  point  source  at  infinity,  we  can  use  the  approximations: 


exp(-ik|r-r' j )  m 
4tt  I  r-r 1  j 


exp(-ik(r,-r.ylr'i'| 
4tiT  ' 


(36) 
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c  r 


±r  expi-ifcc*  -  ^  *)]  =  e’1C? 


cc'  * - •*  kr'  as  seen  in  Appendix  B 

cr '  w 


Finally 


exp( ikr*v'r' )  = 


2  mi  l  ^•(-i'.cosOoJsinf-f c.n) 

x  R^(-ic,U)  cosmU-*’)  (39) 


r-v'r1  =  r  •[sinoQcos<{*'  ex  +  si neQsi '  ey  +  cos90  e2] 
fVr'  =  x  si no0cos<f> 1  +  y  sine0sin<!>‘  +  z  cose0 


For  a  wave  incident  from  the  positive  x  half  x-z  plane,  <$>'=0 
and  we  obtain  the  expansion 

*  00  2-Onm 

exp[i k(x  sin<j0+z  cosoQ)]  =  2  [  l  -m— -  in  S  (-ic.cose  ) 

m=0  n=0  mn 


cos  m$  .  (41) 


We  define  the  coefficients,  Ymn(°0)»  of  the  scalar  wave  functions, 
as  in  Flammer[10],  Equation  (7): 

(2-‘W  . 

fmn' °o  ^  ~  ^  Nmn  ’  Smn(-ic,coso0)  (42) 


<n  r , 

exp[ik(x  sino0  +  z  coso0)]  =  £  £  Yran(e  )*^]  }(n,C,i>) 

m=0  n=m  mn 


2Q 


B-  Vector  Plane  Wave  Expansion 

(29).  ^  liSt  thG  obtained  formu1as>  with  the  notations  of  Equation 


ey  exp[ik(x  sineQ  +  z  cose0)] 


L  1  Ymn(fU  MPx(1)-n,f,o) 


ik  coseo  m=0  n=m  WV  \ 


(-ex  cose0  +  e2  sine0)exp[ik(x  sine0  +  z  coseQ)] 

-j  oo  oo 

Tk  ^  Ymn^eo^ 
m=0  n=m  mn 

ey  exp[ik(x  sineQ  +  z  cose0)] 


ik  sineo  m=0  n In  ^  ){ruK>*) 


(_ex  C0S90  +  ez  sine0)exp[ik(x  sine0  +  z  cose0)] 

**»  ou 

5  '  FE550-  v  1  W8o)  (i.c.o  (4 

K  COS m=g  n=m  Cmn 

A  oo  oo 

eJ  exp[ik(x  sin«0  +  z  cose0)]  -  i  J  J  r„(9„)  if/1’ >(n,c,*) 

m=0  n=m  mn 

(4 

("ex  cos9o  +  ®z  Sin90)exp[ik(x  sine0  +  z  coseQ)] 

"  rrk;  J0  X  Y»"<e°)  <4! 
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C.  Normal  Incidence  Case,  o0=0 

In  this  section,  we  will  study  the  behavior  of  Equations  (44)  to 
(49)  in  the  case  of  a  plane  wave  normally  incident  on  the  top  of  the 
disk.  As  we  have  seen  in  Formula  (15),  the  angular  functions  of  order 
greater  than  zero  are  equal  to  zero  for  unity  argument. 

m£l  Smn(-ic ,cose0)  =  0  for  0O  =  0. 

From  Equation  (42),  we  deduce  that  yQn(0)  is  the  only  non-zero 
expansion  coefficient.  The  expansions  Equations  (44),  (45),  (47)  and 
(48)  can  be  readily  seen  to  be  single  summations  on  n,  0<n<«>  with  m=0. 
The  expansions  Equations  (46)  and  (49),  however,  involve  a  l/sine0 
factor  that  leads  to  indeterminacy  of  the  limit  when  e0  tends  towards 
zero. 

As  we  can  see  from  the  definition  of  Smn(-ic,n),  Equation  (8), 
the  angular  function  is  led  by  a  ( 1 -n2)m/2  factor.  Therefore,  from 
Equation  (42),  we  can  write: 

W9o)  “  0-cos2e0)m/2  =  sinme0  (50) 


W°o>  ■ 


We  can  then  conclude  that  the  value  of  the  limit  when  e0  goes 
to  zero  for  the  different  values  of  m. 

.  .  Ymn^°o^  _  q 
For  m>l  1 lm  “rTnn  u 
6o->0  s  o 

For  m=l ,  we  have  seen  in  Equation  (16)  that: 

11.  !  ,  <r+l)(r«)  d’"Hc, 

0  >0  sin°0  2  r=0,l 
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For  m=0,  YQn(o0)/sine(;)  goes  to  infinity  when  e0  approaches  zero. 

We  will  show  that,  despite  the  fact  that  each  coefficient  diverges,  the 
partial  summations  over  n  for  m=0  in  Equations  (46)  and  (49)  are  globally 
equal  to  zero.  This  will  lead  to  two  dependence  relationships  for  the 
spheroidal  vector  wave  functions.  Consider  first  the  case  of  Equation 
(46).  Because  of  an  m  coefficient,  the  n-  and  ^-component  of  Me5^'  are 
equal  to  zero.  We  prove  that  the  ^-component  of  the  partial  summStion 
for  m=0  in  Equation  (46)  is  also  equal  to  zero;  this  is  expressed  in  the 
following  equation: 


TP  S0n(-ic,n)  R^(-ic,U)-£  S0n(-ic,n)  (52) 


<Lr 
de  K0n 


{1)(-ic,ic)]=0 


'J 


This  must  be  valid  for  any  pair  (n»0-  In  order  to  prove  this 
identity,  we  must  go  back  to  the  scalar  plane  wave,  Equation  (41),  and 
specialize  to  our  normal  incidence  case. 


exp[ikz] 


°°  (11 

I  W°)  W'ic‘n)  Ron  (“ic’ic) 


(53) 


From  Appendix  B,  we  know  that  z  can  be  expressed  very  simply  in 
terms  of  n  and  F/. 

z  =  an?  -  (54) 

Taking  the  derivative  of  Equation  (53)  with  respect  to  n  and  f, 
separately,  we  obtain  two  expansions  as  follows 


~  [exp(ikz)]  =  ikac  exp(ikz)  =  l  Y0n^  S6n^“ic,n^  R0n^"ic,i^ 
an  n=0 

(55) 


Jr  [exp(ikz)]  =  ikan  exp(ikz)  =  [  Y0n(0)  $0nHc,n) 

^  n=0 


(-ic.ic) 

(56) 


Multiplying  Equation  (55)  by  n  and  Equation  (56)  by  -q  and  adding 
them  leads  directly  to  Formula  (52).  We  have  shown  that,  for  any  Dair 
(n.c): 


>0n<°>  Meoi1)<n't'+> 


0 


(57) 
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Taking  the  curl  of  this  equation  leads  to: 


Jo  W01  S4l'!("-f-*)  =  “  (58) 

To  suppress  the  indeterminacy  of  the  limit  of  the  ratio  of  Equations 
(57)  and  (58)  to  sinoQ  when  e0  goes  to  zero,  we  will  apply  L'Hospital's 
rule  to  each  component  of  those  equations: 


lim  -J—  l  y  (e0)[FL 
n  .n  sine„  nEn  On  °/L  B 


o  >0  sine0  n=0 


=  lim  l 


?in  S0n(ic,coseo)(-sine  ) 

-  1  — -  ■  ” 1  -  - - -  r  n 


o0  >0  n=0  "On 


where  [F]e  represents  the  B-component  of  either  Mpz^  or  N  z^ 


B=n»t , 


T5U//1o?  conclude  that,  in  the  normal  incidence  case,  the  expansions 
(46)  and  (49)  can  be  expressed  as  single  summations  on  n,  l<n<»  with  m=l 
[;lammer[7]  gives  a  vectorial  proof  of  Equation  (57)  by  taking  the  curl  of 
ez  exp(ikz).  Equation  (53),  which  is  zero  as  the  curl  of  a  z-directed 
vector  whose  component  is  a  function  of  z  only. 

For  convenience  in  the  study  of  the  normal  incidence  case,  we  de¬ 
fine,  as  F!ammer[10],  the  following  coefficients: 

oo 

-On  =  2  *"■'  <  V  d?n(-fc) 
r=0,l 

■  2  KoI  50n(-ic,l) 


I 

i  Y0n 


b,  =  2  ("*'  n:t  T'  inm  d,n(-1c) 

In  In  r4_,  r!  r  <  '<=) 

The  vector  wave  expansions  are  then,  for  normal  incidence: 


e  exp(ikz)  =  1  f  a  M  x(1 
y  k  n=Q  On  eQn 
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(63) 


In  the  previous  two  sections,  we  have  derived  various  relation¬ 
ships  between  spheroidal  scalar  functions,  and  between  spheroidal  vector 
wave  functions.  Some  limiting  cases  will  allow  us  to  show  some  depend¬ 
ence  relationships,  as  follows. 


Equations  (55)  for  £=0  shows  that  the  derivatives  of  SQ  (-ic,n) 
are  not  independent  and  therefore  not  orthogonal  as  the  original  func¬ 
tions  were. 


Equations  (57)  and  (58)  show  the  lack  of  independence  of  the  func¬ 
tions  M6Qn  ,  and  Neon^’  resPectively.  Equations  (44)  and  (47)  ^multi¬ 
plied  by  cos0o,  lead  to  dependence  relationships  for  MexO),  and  Ngx(  ' 

respectively,  when  approaches  tt/2.  That  kind  of  rela?ionship  cannbe 
derived  for  each  family  of  vector  wave  functions  of  index  i=l.  We  do 
not  deal  with  independent  sets  as  in  the  spherical  case.  In  the  next 
two  parts,  we  will  derive  two  different  solutions  for  the  scattering 
problem.  In  both  cases,  the  fields  will  be  expanded  in  spheroidal  vec¬ 
tor  wave  functions.  For  the  incident  fields,  we  will  use  the  expansions 
derived  in  the  previous  two  sections.  We  will  not  be  able  to  consider 
individual  terms  in  the  summation  expressing  the  incident  field  to  be 
modes  as  we  do  in  simpler  problems,  where  orthogonality  properties  exist. 
To  each  vector  wave  function  of  the  expansion  of  the  incident  field  will 
correspond  not  only  one  vector  wave  function  in  the  scattered  field  but 
a  summation  of  the  scattered  counterparts  of  all  the  vector  wave  func¬ 
tions  appearing  in  the  dependence  relationship  satisfied  by  the  initial 
wave  function.  This  will  be  shown  in  part  VI-C.  We  have  not,  however, 
proved  that  the  spheroidal  vector  wave  function  of  index  i=4  are 


25 


dependent  since  the  relationships  obtained  for  index  1  are  not  appl 
cable  to  functions  of  index  i=4. 

In  the  following  chapter,  we  will  rederive  Flammer's  solution 
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CHAPTER  V 

FLAMMER'S  SOLUTION 


Unlike  Meixner,  F1ammer[7]  deals  directly  with  the  electric  and 
magnetic  fields  in  order  to  find  a  solution  to  the  scattering  problem. 
This  approach  allows  simpler,  more  physically  interpretable  boundary 
conditions,  but  at  the  expense  of  more  complicated  expansions.  The 
major  problem  encountered  in  the  scattering  by  a  disk  is  that  satis¬ 
fying  the  boundary  conditions  on  the  surface  of  the  disk  leads  to  solu¬ 
tions  that  do  not  necessarily  satisfy  the  edge  condition.  The  edge  con¬ 
dition  puts  restrictions  on  the  behavior  of  the  components  of  the  total 
electromagnetic  field  in  the  vicinity  of  the  edge  as  explained  in  Chap¬ 
ter  II  of  this  work.  For  the  electromagnetic  energy  density  to  be  inte¬ 
grate  in  the  vicinity  of  the  edge  all  components  of  the  fields  vary  at 
most  as  s'  '  where  s  is  the  distance  from  the  edge,  while  the  ^-compon¬ 
ent  of  the  total  magnetic  field  remains  finite  and  the  <(> -component  of 
the  total  electric  field  goes  to  zero  as  s''  .  A  very  remarkable  point 
in  Flammer's  solution  is  that  the  satisfaction  of  the  condition  on  the 
(fi-component  of  the  electric  field  directly  ensures  the  satisfaction  of 
the  edge  condition. 

Whatever  the  expansion  of  the  electric  field  of  the  incident  plane 
wave  in  spheroidal  vector  wave  functions  of  index  1,  the  summation  of 
spheroidal  vector  wave  functions  of  index  4  which  satisfies  the  E-field 
boundary  conditions  on  the  surface  of  the  disk  has  a  singularity  of 
order  s"^  in  its  ^-component.  Flammer  chose  to  express  the  incident 
plane  wave  as  sum  of  two  different  expansions.  The  relative  weights  of 
each  expansion  are  then  chosen  in  such  a  way  that  the  coefficients  of 
the  singularity  in  the  ^-component  of  the  corresponding  E-field  expan¬ 
sions  of  index  4  cancel  each  other.  With  the  proper  choice  of  the  two 
plane  wave  expansions,  the  above  procedure  also  leads  to  the  satisfac¬ 
tion  of  the  whole  edge  condition. 

First  we  will  study  the  normal  incidence,  perpendicular  polariza¬ 
tion  bistatic  case  because  of  its  greater  simplicity.  Perpendicular 
polarization  is  here  understood  as  defined  in  Chapter  II,  considering 
the  x-z  plane  as  plane  of  incidence.  The  problems  encountered  in  ex¬ 
tending  the  method  to  the  arbitrary  incidence,  arbitrary  polarization 
bistatic  case  will  be  noted,  and  the  results  of  a  numerical  test  of  the 
far-field  for  the  normal  incidence  case  will  be  presented. 
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A. 


Normal  Incidence,  Perpendicular 
Polarization  Bi static  Case 


This^  corresponds  to  the  case  ot=n/2,  E1|)ey,  and  FP||e„.  The 
incident  E-field  can  be  expanded  according  to  Equations  (62;,  (64)  or 
(66)  and  the  corresponding  H-field  according  to  Equations  (65),  (67) 
or  (63),  respectively.  To  represent  the  outgoing  scattered  waves,  we 
use  vector  wave  functions  of  index  4  which  match  spherical  waves  when 
i  becomes  large. 

We  will  successively  consider  the  boundary  conditions  on  the  sur¬ 
face  of  the  disk,  an  expansion  of  the  ^-component  of  the  electric  field 
in  the  vicinity  of  the  edge,  the  [P],  condition,  and  the  final  solution 
where  we  will  check  the  satisfaction^  the  edge  condition. 

1 .  Field  boundary  conditions 
on  the  surface  of  the  disk 


We  will  successively  consider  the  three  possible  expansions  and 
the  scattered  fields  they  lead  to.  Despite  the  lack  of  orthogonality 
of  the  spheroidal  vector  wave  functions  of  index  1,  we  will  use  a  term 
by  term,  function  by  function,  matching  technique  to  obtain  the  expan¬ 
sion  of  spheroidal  vector  wave  functions  of  order  4  which  satisfies 
the  boundary  condition  on  the  surface  of  the  disk  for  given  incident 
waveexpansion.  These  expansions  will  be  called  "reflected"  fields  P" 
and  Hr.  The  reflected  fields  are  not  equal_.to  the  scattered  fields 
unless  the  edge  condition  is  satisfied  by  E"'+E'r  and  H1+Hr,  but  they 
satisfy  the  boundary  conditions  on  the  surface  of  the  disk.  The  bound 
ary  condition  is  given  in  Equation  (4): 

[f+E5],  =  [f+E5]^  =  [hVh5^  =  0 


The  term  by  term  matching  technique  does  not  lead  directly  to  the  scat¬ 
tered  field  since  the  spheroidal  vector  wave  functions  of  index  1  are 
not  independent.  We  now  express  the  reflected  fields  for  the  three 
possible  expansions  of  the  incident  fields. 


-  x(l ) 

Mq  v  '  expansion 
_ 

From  Equations  (62)  and  (65), 
E 


n  .  2  v  ,  n-  xO) 
'  k  n=0  0n  e0n 


rji  _  ,•  H°  V  a  N  x  ^ 1  ^ 

H  -  '  r  „i0  a0n  %n 


where  H0=E0/Z0  and  Z0  is  the  free  space  impedance. 


(68) 


i 

t 
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From  Equations  (30)  to  (33),  and  Equation  (22),  we  see  that,  for 

n  even,  the  components  of  M  *0)  tangential  to  the  disk  and  the  normal 
_  F7i  \  eQn 

component  of  Ne^  '  are  equal  to  zero.  Therefore,  in  Equation  (68), 

the  vector  wave  functions  for  n  even  do  not  lead  to  any  reflected. field. 
For  n  odd,  the  above  three  components  are  all  proportional  to  RX^(-ic, 
iO),  and  the  same  components  of  the  functions  of  index  4  are  propor¬ 
tional  to  R^J)(-ic,iO).  Since  the  only  difference  between  a  spheroidal 
vector  wave  function  of  index  1  and  its  counterpart  of  index  4  is  their 
radial  dependence,  it  is  possible  to  cancel  a  component  of  one  by  that 
of  the  other  on  the  surface  of  the  disk  with  appropriate  coefficients 
in  order  to  satisfy  the  boundary  condition. 


F(])(n,CM) 


R0n  (-ic>i0) 
Rg^  ( -  i  c ,  i  0 ) 


F^4)  (n,0,<!>)  =  0  for  0<|n|lO 


where 


F^1 )(  n  £  <i>)  =  I'M  h  FM  h  TN  x(i)l 

i-  ln.5.W  *-Me0n  Jn  »  L%n  »  LINe0n  h 


The  term  by  term  cancellation  method  leads  to  the  following 
reflected  fields  which  satisfy  the  boundary  conditions  on  the  surface 
of  the  disk 


FT  _  Si  y  .  R0n1'(-’c-,0)  jT*<4) 

k  n-1  °"  R‘JI'(-ic,iO)  e°" 


( -ic,iO)  _ 


nr  = _ 0  y i  On _ _  «7 

k  n=l  0n  R^'f-iciO)  e°" 


where  the  prime  over  the  summation  means  that  the  index  of  the  summation 
varies  by  increments  of  2  from  its  initial  value.  This  notation  will  be 
used  throughout  the  rest  of  this  work.  In  the  same  way,  for  the  other 
two  expansions,  we  obtain  the  following  formulations. 


»*(!> 


expansion 


From  Equations  (64)  and  (67), 
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(72) 


Frotn  Equations  (30)  to  (33),  and  Equation  (22),  we  see  that  for  n-1 
odd,  the  components  of  Mg^  ^  tangential  to  the  disk  and  the  normal 

component  of  are  equal  to  zero.  For  n-1  even,  the  same  com¬ 

ponents  are  proportional  to  Rj^(-ic,iO)  and  the  corresponding  compon¬ 
ents  of  the  functions  of  index  4  to  r(^ ( -ic, iO) . 

I  n 

Using  the  same  term  by  term  matching  technique,  we  obtain  the 
following  reflected  fields: 

E.  Rp^(-ic,iO)  _  7(4) 

7=r  Ori,  In  \a  '  i  7") ' 

E  -  ■  r  bl"  e,n 


i  H  '» 

TT  I;  "in 
n=l 


R^n)(~ic»iO)  _  z(4) 
-m -  N 


R|n  ("iC’iO) 


—  v(  1 ) 

No  expansion 
K0n _ 

From  Equations  (66)  and  (63), 
fi  =  1^1  y  a  n  h“ 

k  n=o  °n  * 


^  V  a  My(1) 

<  Jna°"%n 


From  Equations  (30)  and  (33),  and  Equation  (22),  we  can  see  that  for 

n  odd,  the  components  of  N  ^  tangential  to  the  disk  and  the  normal 
_  v(i l  e0n 

component  of  are  equal  to  zero.  For  n  even,  the  same  compon- 

eOn 

ents  are  proportional  to  Rq^  '  (-ic.. '0y  and  the  corresoonding  components 
of  the  functions  of  index  4  a>-e  proportional  to  Rg^4'(-ic,i0). 

Using  the  term  by  term  matching  technique,  we  obtain: 


ET  =  - 


-°-  P  a  RQn)H<:-i°)  jj  y<«) 
k  n=0  0n  R^J^-ic.iO)  e°n 
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H 


(77) 


V  ?.  a  ^’hc.iq) 

k  n=0  0n  R^I-ic.iO) 


m  y(4) 

e0n 


All  three  pairs  of  incident  and  reflected  fields  satisfy  the 
boundary  condition  on  the  disk,  excluding  the  edge.  We  will  study  the 
behavior  of  [F’+F1'’]^  at  the  edge  in  the  following  section  for  each  case. 


2.  [E]^  condition  at  the  edge 

We  will  first  give  the  expansion  of  [F’+Fr].  in  the  vicinity  of 
the  edge,  and  then  present  Flammer's  method. 


a.  Behavior  of  [F1+Er]i 

_ _ _ i 

at  the  edge 

The  edge  is  located  at  the  coordinates  n=0  and  £=0.  We  will  use 
the  small  argument  expansions  for  the. angular  and  radial  functions  to 
obtain  a  power  series  expansion  of  [Ei+Er](j).  In  each  summation,  we 
will  only  retain  the  diverging  and  constant  terms,  when  they  exist,  and 
the  converging  term  of  smallest  order  at  the  edge.  We  will  use  the 
following  expansions: 


(l-o2)1/2  =1  .  1  n2  +  O(n^)  for  n  and  £  small 
d-n2)-l/2  =  t  +  ^n2  +  0(n4)  Jn|«  1 

(1h2)1/2  =  1  +  1.  c2  +  0(54)  5  «  i 

(1+£2)'1/2  =  1  -  1  C2  +  0(£4) 


We  will  now  consider  the  three  cases. 


—  x(l ) 

Case  of  the  Me  expansion 


(78) 


A 


For  this  case  only  we  will  give  the  whole  procedure: 


t 
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[£l+E,%  ;  k  )•'  a0n  d  2  2  ( 1  -  n2 )  n  ( -  a  Qn )  ,1  +  AQn  j- ^ 

n=0  +n  L  a 


+ 


2“  On  / 


.  I  y.  a  2  _J _ 

k  i  aOn  d  2  2 

K  n=l  u  f3  +n 


+  n(l+#:2)n  U  +  §-  (;0n 


J-  v 1  ,  -  un 

k  Z  aOn  d  -T4 
n=l 


^  n(l+C2)f.  xQnj  Rjp(-ic,iO)S0n(-ic,0)cos-, 
[f.O-n2)(l  -f  Uon-2)? 

+  §“  R&^'(-ic.10)S|jn(-ic,0)cos* 

* 

R^}  ‘  (-Tc,iO)  ,  r  /  ,2  \ 

HKJ 


1  v, 2  RyMcjo)  1  rn  2 


R0n  ^'ic>i°)  ?2+r|2 


‘  2  ^A0n~2)) 


+  nO+fr)n£  A, 


i  S'  ( - i c , 0 ) 


c  RQn'  (-ic.iO) 


+  S(l-n2) 


(|  -  ^^On-2)^ 


T+oO+S2)  +|_(>^n_2) 


Rq^  ( - i c , i 0 )  x  S'  (-ic,0)  l  cos} 


We  can  see  directly  that  the  n=odd  terms  of  the  incident  wave  are 
cancelled  by  part  of  the  terms  of  the  "reflected"  field  in  this  expan¬ 
sion.  Finally  Equation  (79)  reduces  to: 

*  h  "E  L 

n-u  i 


,V  »=1  c 


l2°-l  fi  *  V  <„V>] 


R0n>,(-iC*i0)  | 

x  cos<f>  (80) 


The  singularity  at  the  edge  comes  from  the  second  term, 
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r  r,  1  "'On 
n=l  C  a°n  ^ 


Snn(-ic.O) 


ROn  ('ic»iO) 


which  does  not  converge  when  c->0  and  n->-0. 

Case  of  the  Mp^1)  expansion 
_ Ijj! _ _ _ 

The  same  procedure,  taking  into  account  the  use  of  the  approxima¬ 
tions  of  Equation  (78)  leads  to  the  following  result. 


■  -  h 


Aln'3 


n  €  l  b1n~-S;n(-ic,0)R)^  (-ic.iO) 
n=2 


00  hn  T  O  O  1  Si  n(“ic>0) 

+  -j-j  I*i  -ID-  1  +  -If-  (nV)  4)  ' - 

nV  n=l  c  L  2  J  Rin  (-1c-i 


0) 


COSO 


(81) 


We  notice  that  the  singularities  in  Equations  (88)  and  (81)  are 
of  the  same  order  C/(n2+C2).  This  is  the  reason  why  Flammer  is  able  to 
cancel  the  singularity  by  taking  the  incident  field  as  a  linear  combi¬ 
nation  of  Equations  (68)  and  (72)  with  proper  relative  weights  of  the 
two  expansions. 


Case  of  the 

_ _ ^On 


expansion 


The  [NeJ(1)] 

angular  and  raOial 
compute.  We  will 


component  involves  second  order  derivatives  of  the 

functions  and  are,  therefore,  more  complicated  to 
use  the  following  expression  for  that  component: 


[n  y(i)i 

e0n  « 


*  '  ka2  ^  (s°n>  R^") 

+  (,-n2>  ^  <V  +  «  s0n  h  (Rto> 


+  <'*c2)s0n 


(82) 


Using  Equations  (13),  (14),  (23),  (24),  (26)  and  (27),  we  finally 
obtain  for  the  tangential  electric  field  at  the  rim: 
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[F+r'i 


k2d2 


cost 


:  3( ‘0n_2) 

+  lrr'  ^  a0n - 2 -  s0n(-ic>°) 

n=  I 


(1)' 

Ron  ( ” i c  > i 0 )  + 


(A0n-2) 


n2+f2  n=0  - 
*n„  S0n(-ic,0) 


i1  +  ^(-o„^2)  - 


an 

)  On 
c 


«On  (-ic^0) 


(83) 


These  three  expansions  of  the  ^-component  of  E1+rr  for  different 
incident  plane  wave  expansions  behave  in  the  same  way  in  the  vicinity 
of  the  edge.  They  each  have  a  singularity  which  comes  from  the  ill-de¬ 
fined  spheroidal  vector  wave  functions  at  the  edge.  We  have  seen  by 
comparing  the  power  series  expansions  of  r(  j/(-ic,i;)  and  its  derivatives 

for  small  argument  ,  Equations  (23)  and  (24),  to  those  of  (-ic,i;) , 
Equations  (26)  and  (27),  that  part  of  the  expansions  of  the  ?0nctions  of 
index  4  is  identical  to  the  whole  expansion  of  the  function  of  index  1 
except  for  the  change  of  index.  The  term  by  term  matching  technique 
exactly  cancels  these  terms  between  the  incident  and  reflected  fields. 
Therefore,  in  Equations  (80),  (81)  and  (83),  the  second  summation  con¬ 
taining  the  singularity  comes  from  the  terms  in  the  expansions  of  the 
radial  functions  of  index  4  that  do  not  disappear  in  the  term  bv  term 
matching.  Their  coefficient  r/( n2+f2)  shows  that,  as  expected,  E1+Er  is 
equal  to  zero  on  the  surface  of  the  disk.  Furthermore,  the  first  sum¬ 
mation  in  Equations  (80),  (81)  and  (83)  corresponds  to  the  part  of  the 
incident  plane  wave  expansion  that  does  not  lead  to  a  reflected  field 
since  it  is  equal  to  zero  on  the  surface  of  the  disk.  This  summation  is 
also  equal  to  zero  at  the  rim  as  the  factor  tv;  shows. 


b.  Flammer's  method  . 

In  order  to  satisfy  the  [E] .  condition,  Flammer[7]  simply  chose 
to  express  the  incident  wave  as  the  sum  of  two  different  expansions 
whose  weights  would  be  adjusted  so  that  the  singularity  in  the  compon¬ 
ent  of  the  electric  field  at  the  edge  would  disappear.  In  order  to 
satisfy  the  whole  edge  condition  by  the  procedure,  it  is  necessary  to 
choose  the  two  expansions  so  that  the  other  components  of  the  fields 
haye  the  right  behavior  near  the  edge.  Tie  n-  and  ^-components  of 

^  have  higher  order  singulari ties  than  those  of  Mp*'  '  and  Mp*^, 
e0n  3  3  e0n  eln 

while  the  latter  have  the  same  singularities,  and  so  do  their  magnetic 
field  counterparts  N^1 )  and  N^1).  We  will  show  that  the  choice  of 
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the  two  expansions  in  terms  of  M  Equation  (68),  and 

e0n  eln 

Equation  (72),  leads  to  a  satisfactory  solution. 


1  a°n<1)+*z  l  bln^1} 


with  3X+3Z=1  in  order  to  have  |F|=E0.  '* 

The  scattered  field  will  be: 

n  1  • 

^  Eo  „  R0n  Hc.iO)  jt  x(4) 

E  "  k  P  7  a°"  R^)'(-ic.iO)  e°" 

°°  Rln^("1c»iO)  -  /  \ 

+  h  V  bln  “TO - 7  Mei(4)  ■  (85) 

n=l  Rln  (“iC’iO)  ln 

From  Equations  (80)  and  (81),  we  see  that  the  coefficient  of  the 
singularity  of  the  expansion  of  [F+Es]^  is: 

o  “  i  (-ic,0)  „  v.  1  .  sin(_1c>0^ 

R  y  la„  0"  .... _ —  +  Bz  Z  r  b-i  -nn -  coS4> 

kd  Bx  A  c  a0n  7(4)  .  TO  n=l  c  1n  r|4 


R0n  (-ic.iO) 


n  ( -  i  c ,  i  0 ) 


In  order  to  set  this  coefficient  equal  to  zero,  we  must  choose 
8X  and  3Z  as  follows: 


With  ex  and  8  defined  as  above,  F+F  satisfies  the  [EF  condi¬ 
tion  at  the  edge.  Equation  (85)  gives,  then,  the  bistatic  scattered 
field  of. the  disk  with  normal  incidence  and  perpendicular  polarization 
(i.e.,  F||£y).  In  the  next  section,  we  will  verify  that  Es  satisfies 
the  whole  edge  condition,  underlining  the  convenience  of  Flammer's  met¬ 
hod  which  satisfies  the  edge  condition  by  dealing  only  with  the  condi¬ 
tion  on  the  <)> -component  of  the  electric  field. 
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3.  The  bistatic  scattered  field 


The  total  scattered  electromagnetic  field  obtained  with  Flammer's 
solution  is,  in  the  normal  incidence  case: 


= 


ft5 


\*i  i: 


n=l 


(’ll1 

R0n  (_ic,iO)  -  x(4) 

a°"  i^Hcao) 


x  -  Rp^-ic.iO) 

+  [iz  l'  bln  ~m - Mei 

n=l  Rln  (”ic*iO)  ln 


iH, 


<  i: * 

n=l 


(l)1 

R0n  Hc.iO)  x(4) 


On  JJT 


Ron' 


e0n 


+  y.  b  R]n>(-i':-1°>  jyzd) 
z  n=l  ln  R[^(-ic,iO)  6ln 


where  6*  and  6*  are  defined  by  Equation  (87). 


(88) 


First  we  will  verify  that  the  above  fields  satisfy  the  whole  edge 
condition  in  the  vicinity  of  the  edge.  The  behavior  of  each  component 
of  the  fields  will  be  expressed  in  terms  of  local  coordinates  that  have 
been  defined  in  Appendix  A.  Second,  the  bistatic  scattered  far  field 
will  be  computed  in  spherical  coordinates. 


a.  Total  field  in  the  vicinity  of  the  edge 

The  edge  condition  is  expressed  in  terms  of  the  distance  from  the 
edge  to  the  observation  point.  We  will  introduce  the  local  coordinate 
system,  valid  for  n  and  f  small,  shown  in  Figure  5-1.  We  use  the  follow¬ 
ing  notation: 

s  =  distance  from  the  edge,  s<<a 

t  =  angle  between  the  top  surface  of  the  disk,  n>0, 
and  the  direction  from  the  edge  to  the  observation 
point,  0<t<2n. 

■I-  =  same  as  in  spheroidal  coordinates. 
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Figure  5-1.  (s,t,<ji)  coordinate  system. 


The  correspondance  between  spheroidal  coordinates  and  these  co¬ 
ordinates  is  derived  in  Appendix  A. 

's  =  j  (n2+C2) 

4 

2  2 

cost  =  \—y  (89) 

n  H 

>s,nt  =  -fe 

n  +£ 

or,  conversely. 


sin(t/2) 


(90) 


The  disk  is  characterized  by  t=0. 

In  order  to  obtain  the  components  of  the  fields  in  the  new  coordi 
nate  system,  we  first  expand  the  radial  and  angular  functions  using 
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their  small  argument  approximations  and  then  transform  the  terms  of  lower 
order  with  the  above  formulas. 

For  the  n-component  of  E=E1+Er,  the  singularity  in  [Me*^]n  and 
[Me^^]n  is  ( o 2+f.2)"1/2,  or  s-l/2t  ancj  has  a  satisfactory  behavior 

near  the  edge.  More  precisely,  we  have 

=  E0  ej  1  r  a0„  (-f)  sint/2  s0n(-^,0)RO>,,(-ic.i0)s1n« 

K  n=0  x  / 
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x  r(^)  (-ic,iO)sini}> 
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.  ,  "■  2  IT  r  R{l*(-ic,iO) 

E°  62  R  Ji  S’"'-iC>0)  [  ^7777] 


x  (-ic,iO)sin<|> 
In 


The  first  summation  is  equal  to  zero  on  the  surface  of  the  disk 
but  apparently  only  there.  However,  if  we  differentiate  Equation  (56) 
again  by  £,  we  obtain 


C  ^ 

— o  [exp(ikz)>  -k2a2n2  exp(ikz)  =  l  yn  (0)  S0n(-ic,r,) 

3f/  n=0  un 

' 1  \« 

x  R0n  HC>U) 


At  the  edge,  where  f,=0  and  n=0.  Equation  (92)  becomes: 


>:•  >0n(0)  S0n<-ic-°>R^,"(-fc.f0)  --  0 
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The  terms  of  this  summation  for  n  odd  disappear  since,  for  n  odd, 
SQn(-ic,0)  is  equal  to  zero,  Equation  (14).  Equation  (93)  is  the  coef¬ 
ficient  of  sint/2  in  the  first  summation  of  Equation  (91),  recalling 
that  by  definition  a0n=-iA0n(O) •  Final ly  we  see  that  Equation  (91) 

contains  no  term  larger  than  s^/2  and  therefore,  we  have  at  the  edge: 
[E]n  =  0(5^2)  sin<f>  sint  (94) 

For  the  e-component  of  E,  we  see  that  the  singularity  of 

and  [M  *(^)]  -js  due  to  a  factor  (l/(n^+C^) ^) »  or  s"^2>  uee  the 
definitiBn  of  the  functions  in  Flammer[10])  and  therefore  [FL  satisfy 
the  edge  condition  for  this  component. 


[E]  =  E0  ^  I1  a0n  I  cost/2  S0n ( - i c , 0)  Rjn'(-ic,iO)sin* 
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+  E0  IT  K  bln  %72  ^  sin  t/2  ~TTT, - T  sin^  (95) 

k  n=l  cdJ/2  R) ^; ( - i c , i 0 ) 

The  second  term  varies  as  (s)"^2  anc|  -js  therefore  predominant. 


n=l  cd' 


[E]r  =  0( s~  ‘/2)  s-j, 


[EL  has  been  computed  to  match  the  corresponding  condition.  We 
obtain  from  Equations  (80)  and  (81): 
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Thus,  as  expected. 


[E]*  =  Ofs1/2)  cos<f>  sint/2 


(98) 


The  total  electric  field  obtained  from  Equations  (84)  and  (85) 
satisfies  the  edge  condition.  We  now  need  to  check  the  components  of 
the  total  magnetic  field,  which  are  expressed  in  terms  of  the  functions 
and  .  Because  of  the  greater  complexity  of  these  functions, 

we  will  only  determine  their  order  of  variation.  Since  the  order  of  the 
singularities  of  each  component  is  not  apparent  directly  from  the  defini¬ 
tion  of  the  functions,  we  will  write  out  the  diverging  terms  for  each 
component  as  function  of  n  and  £,  and  then  we  can  write  the  components 
of  the  total  magnetic  field  H.  Since  we  will  only  be  looking  at  singu¬ 
larities,  we  will  not  include  in  the  following  expressions,  the  summa¬ 
tions  corresponding  to  the  incident  waves  since  they  are  well  defined 
at  the  edge,  that  is  to  say,  of  order  1  or  smaller. 

For  the  n-component  of  H,  we  have  the  following  singularities  at 
the  edge: 


For  n  odd. 
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From  Equations  (89)  and  (90)  we  find  that: 
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Again  by  definition  of  ex  and  ez,  the  coefficient  of  the  improper 
singularity  in  [H]f,  normal  to  the  disk,  disappears.  Thus, 


[Hi  =  0(s“1/2)  cost,  sint 


(106) 


For  the  ^-component  of  H.  the  singularities  of  the  vector  wave 
function  involved  are  of  order  s"^2.  They  are,  however,  cancelled 
by  the  coefficients  introduced  to  satisfy  the  edge  condition  on  [E]^. 

For  n  odd, 


Ck<j2  R0n^  ( ~ 1  c , 1 0 )  n2+?2 


(107) 
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ln  *  ckd2  Rj^(-ic.iO)  n  +r 


where 


n  =  1  17 


cost/2 


From  the  definition  of  the  scattered  magnetic  field  and  from 
Equation  (107),  we  obtain  for  the  total  field: 
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sin<f>+0(  1 ) 
(108) 


By  definition  of  B*  an<l  6z»  the  coefficient  of  s'^/2  disappears 
and  finally 


[HL  =  0(1)  sin* 


(109) 


Therefore,  the  total  magnetic  field  H  satisfies  the  edge  condition. 
We  have  thus  verified  that  the  scattered  field  obtained  by  Flammer  by 
matching  the  4>-component  of  the  electric  field  at  the  rim  satisfies  the 
whole  edge  condition. 

We  summarize  the  different  results: 


[E]n  =  0(s1/2),  [E]c  =  0(s'1/2),  [E]*  =  0(s1/2) 
[H]r|  =  0(s-1/2)>  [H]c  =  0(s-1/2),  [TT] (j,  =  0(1) 


The  formal  solution  for  the  bistatic,  normal  incidence,  perpendi 
cular  polarization  case  satisfies  the  boundary  conditions  on  the  disk 
including  the  edge  and  is  therefore  proved.  In  the  next  section  we 
will  compute  the  scattered  far-field  for  the  above  case. 


b.  Bistatic  scattered  far-field 
in  normal  incidence 


The  bistatic  scattered  electric  far-field  for  normal  incidence 
and  perpendicular  polarization  can  be  obtained  from  Equation  (88)  by 
replacing  the  vector  wave  functions  by  their  large  radial  argument 
approximations,  listed  in  Equation  (34).  From  Appendix  A,  we  will  use 
the  approximation  of  the  coordinate  system  for  large  £. 
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Thus  the  scattered  field  becomes: 
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(112) 
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We  can  easily  check  that  E^H3  is  parallel  to  er  since  we  are  in 
the  far  field  region. 

The  electric  field  given  by  Equation  (112)  has  been  programmed, 
and  the  results  are  compared  with  those  obtained  from  the  program  de¬ 
rived  by  Hodge[9]  with  Meixner's  solution  in  a  later  section.  We  will 
now  study  FI  aimer's  approach  to  the  arbitrary  incidence,  arbitrary  polar¬ 
ization  bistatic  case. 


B.  Arbitrary  Incidence 

In  the  previous  section  the  disk  scattering  for  a  normally  inci¬ 
dent  plane  wave  with  perpendicular  polarization  has  been  solved.  The 
general  solution  for  arbitrary  incidence  can  be  attempted  following  the 
same  approach. 

The  general  formula  for  the  incident  field,  Equation  (3),  can  be 
separated  into  two  components 


E^  =  E^(cosa  E1'11  +  sina  E11) 

where 

E1"  =  (coseo  ax  -  sin8Q  az)e-1^  ‘r 
E1^  =  ay  e-iRl*r 


(114) 


Let  us  define  E3!  as  the  scattered  field  corresponding  to  the  in¬ 
cident  field  EU.  The  general  scattered  field  will  be  obtained  from 
and  as  follows: 


E3  =  E0(cos«  E3"  +  sina  ESi) 


(115) 


In. order  tp  follow  the  method  of  the  previous  sections,  we  will 
expand  E1"  and  E11  as  sums  of, two  different  expansions.  We  will  Uise 
Equations  (44)  and  (46)  for  E1-1,  and  Equations  (47)  and  (49)  for  E111. 
The  H-fields  can  be  expressed  with  the  corresponding  spheroidal  wave 
functions. 
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where  +  u*  =  1 ,  +  6^  =  1 • 

Unlike  the  case  of  normal  incidence  these  expansions  involve  all 
the  values  of  m  and,  therefore,  all  the  different  ^.-dependences .  The 
reflected  fields  are  calculated  using  the  same  term  by  term  matching 
technique  as  in  the  normal  incidence  case.  From  Equations  (30)  to  (33), 
we  obtain  the  following  expansions. 
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II  II 

For  arbitrary  values  of  and  b£,  the  edge  conditions  is  not  sat¬ 
isfied.  As  in  the  normal  incidence  case,_we  will  try  to  satisfy  the 
edge  condition  by  simply  satisfying  the  [E]^  condition  at  the  edge.  We 
will  study  the  perpendicular  polarization  case  in  more  detail.  The 
method  can  be  exactly  paralleled  for  the  parallel  polarization  case. 

We  will  now  expand  the  4-component  of  the  total  electric  field  in  the 
vicinity  of  the  edge  for  the  perpendicular  polarization  case.  We  need 

the  expansion  of  the  summations  on  and  TTe^^ 1 ) 

In  the  same  way  that  we  derived  the  expansions  of  [E1+Ey"]a  in  the 
vicinity  of  the  edge  for  the  normal  incidence,  perpendicular  polariza- 

_  X  ( 1  ) 

tion  case  for  plane  wave  expansions  in  terms  of  Mp„  ,  Equation  (80), 

—  _ / i \  cun 

and  MefT  ' ,  Equation  (81),  we  obtain  here  in  the  arbitrary  incidence 
case:  ,n 
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(120) 


where  6mn  is  the  standard  Kronecker  symbol. 


47 


x  The  condition  on  [E]+t  ].  is  thdt  it  be  zero  on  the  rim.  and 
.  z  must  then  cancel  the  coefficient  of  the  singularity  .%/(-  ^+T^)in 
Equation  (120).  If  we  reorder  the  coefficient  of  the  singularity  and 
equate  it  to  zero,  we  obtain  the  following  equation: 
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Note  that  the  coefficients  and  have  been  taken  Inside  of  the 
summation  of  index  m. 

Because  of  the  .^-orthogonality,  the  coefficient  of  each  factor 
cosm^  must  be  equal  to  zero.  Because  of  the  relation  3x+:'z=^  >  we  obtain 
for  i',^  and  the  various  following  relationships  that  must  be  satisfied 
simultaneously. 
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for  m=l 
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In  the  limit  when  eQ  tends  to  zero,  these  coefficients  lead  to 
those  of  the  normal  incidence  case. 

For  m>l 
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In  order  to  satisfy  the  [E] ,  condition  at  the  edge,  8*  and  et 

<}>  X  Z 

must  satisfy  Equations  (122),  (123)  and  (124)  for  every  value  of  m. 
According  to  our  definition  in  Equation  (116),  6^  and  8*  are  independ¬ 
ent  of  m  since  they  represent  the  relative  weights  of  the  two  differ¬ 
ent  incident  plane  wave  expansions.  To  see  if  a  pair  (b^S*)  can 
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satisfy  the  [Ej^  condition  at  the  edge,  the  different  expressions  of 

i^,  for  example,  in  Equations  (122)  to  (124)  should  be  computed  for 

various  values  of  m.  However,  since  the  vector  wave  functions  of  index 
i=l  are  not  independent,  the  fact  that  the  coefficients  should  take 
different  values  for  each  m  is  not  a  proof  that  the  solution  is  invalid. 
If  this  happens,  another  approach  introduced  by  Flammer[7]  can  be  used. 

Instead  of  one  pair  (e^.e1),  Flammer  uses  a  whole  set  of  pairs  ( ,e*  ) 

i  1  x  z  xm  zm 

where  B  and  ft  satisfy  the  Equations  (122)  to  (124)  corresponding  to 
xm  zm  , 

the  value  of  m.  Note  that  while  ft,  will  then  be  the  coefficient  of 
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in  the  incident  field,  sv  is  not  the 


weight  of  the  summation  over  n  of  M 


m  xO  ) 


but  of  other  functions  whose 


^.-dependence  is  due  to  one  trinogometric  function  only  per  component, 

—  x  ( "M 

instead  of  product  of  trigonometric  functions  in  the  case  of  Mg'  . 

This  is  due  to  the  reordering  necessary  to  obtain  Equation  (121).  The 

new  functions  Me^ijp  and  Me^in  are  defined  in  Flammer's  book[10]. 

Their  m-index,  m+1  and  m-1  respectively,  represent  the  order  of  their 
-dependence.  In  order  to  write  the  incident  field  used  by  Flammer,  we 
introduce  the  relation 


for  m>l 
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The  incident  field  is  then  equal  to  the  following  formula: 
► 

E  f ,  «>  Y  ( 0  )  /  -i  \ 

lk  m=0  0,11  xm_1  n=m  cos0o  ^  >" 


+  (1  +<5  )  B1  y  Ymn^°°^  M  +0) 

(  60m  Bxm+1  n^m.  cosOQ  em+l,n 


Y  (u  ) 

1  mr\  '  r\  ' 


.  y 

zm  ’  siw> 
m  n=m  o 


mn  o  M-  z ( 1 ) 


(126) 


where  B1  and  [>}  satisfy  the  corresponding  condition.  Equations  (122) 
to  (124)”.  ” 
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A  term  by  term  matching  technique  leads  to  a  reflected  field  which 
satisfies  the  boundary  conditions  on  the  surface  of  the  disk  and  the 

[E],  condition  at  the  edge.  But,  when  these  coefficients  are  utilized 
4 

in  Equation  (126),  it  is  not  clear  whether  these  expansions  still  repre¬ 
sent  plane  waves.  This  question  remains  to  be  answered. 

By  introducing  different  coefficients  for  each  value  of  m,  we  have 
modified  the  initial  plane  wave  expansions.  Since  the  vector  wave  func¬ 
tions  of  index  1  are  not  independent,  the  incident  field  of  Equation 
(126)  may  still  be  a  plane  wave.  This  needs  to  be  checked  for  every 
value  of  e0.  Computations  of  Equation  (126)  on  the  surface  of  the  disk 
should  be  made  to  verify  whether  or  not  Flammer's  solution  is  valid  for 
arbitrary  incidence.  In  his  paper[7],  Flammer  does  not  mention  the 
problem.  It  is  in  no  way  proven  and  a  computational  verification, 
though  not  a  proof,  should  be  made  in  future  work. 


C .  Numerical  Test  of  the  Normal 
Incidence  Case 


A  computer  program  based  on  Equation  (112)  computes  the  scattered 
electric  far  field  in  the  normal  incidence,  perpendicular  polarization, 
bistatic  case.  We  will  successively  introduce  the  variables  computed 
by  the  program,  present  the  computed  data,  and  discuss  the  results. 

1 .  Computer  program 

The  program,  presented  in  Appendix  D,  computes  the  value  of  E, 
where  E  is  defined  by: 


fS  _  Eo  exp(-ikr)  p 
^  "  k  kr  L 


(127) 


The  functions  and  subroutines  presented  in  Reference[9]  were  used 
for  this  purpose.  The  output  data  are,  as  in  Hodge's  program  for 
Meixner's  solution[9]: 

the  magnitude  and  phase  of  the  e-  and  <j>-  components  of  E, 


the  normalized  crpss  sections  of  the  disk  corresponding  to  the 
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e-  and  ^-components  of  ^  as  follows: 


°a^0s,<!)s^  I  4irrC 
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where  a  =  0,41. 


?7|[e]«|2 


(128) 


The  normalization  is  obtained  by  dividing  the  scattering  cross 
section  by  the  surface  area  of  the  disk.  The  magnitude  of  the  cross 
section  in  the  scattering  direction  defined  by  (e  ,<j>  )  is  equal  to 

K2<V*s>+%<V*s»1/2- 

2.  Data 


The  results  obtained  from  Equation  (112)  will  be  compared  to 
computations  made  with  Meixner's  solution.  Meixner's  far-field  solu¬ 
tion  has  been  studied  extensively  by  Hodge[9],  who  wrote  a  computer 
program  which  calculates  the  electric  far-field  for  arbitrary  incidence 
and  polarization.  The  results  of  Hodge's  program  show  very  favorable 
agreement  with  measurements,  with  calculations  based  on  the  small  disk 
approximation,  and  with  Geometrical  Theory  of  Diffraction  (GTD)  compu¬ 
tations.  We  will  accept  that  program  as  a  comparison  source  for  the 
results  obtained  from  Flammer's  solution. 

Various  plots  of  both  solutions  have  been  made  and  we  will  use 
the  following  notation  for  the  curves,  except  where  otherwise  stated: 

- continuous  line  for  Flammer's  solution 

-  dashed  line  for  Meixner's  solution. 

In  the  cases  where  $  =0° ,  H-plane,  or  *=90°,  E-plane,  the  nor¬ 
malized  crqss  section  and  the  phase  will  be  tnat  of  the  non-zero  com¬ 
ponent  of  E,  Oi  when  <t>s~0°  or  when  <j>s=90° ,  respectively.  For  other 
values  of  <j>s,  the  cross  sections  corresponding  to  the  two  components  of 
E  will  be  plotted. 


3.  Resul ts 

From  the  plots  drawn,  we  can  see  generally  good  agreement  between 
both  theories.  We  notice  the  following  tendencies: 

The  concordance  between  the  two  solutions  is  very  good  for  small 
disks  -  ka<2  -.  In  Figures  5.2  to  5.4,  the  low  frequency  ends  of  the 
curves  match  very  well.  Figures  5.5,  5.8,  5.11  and  5.13  show  a  very 
good  agreement  for  both  magnitude  and  phase  in  the  case  ka=2.  For  ka=4. 
Figures  5.6  and  5.14  show  a  reasonably  good_agreement  for  the  E-plane 
cross  section  and  for  the  phase  of  E  in  the  H-plane,  but  the  values  of 
the  H-plane  cross  sections  differ  greatly  as  0]  approaches  90°,  Figure 
5.9.  ka=2  seems  a  reasonable  limit  for  the  low  frequency  region  where 
the  two  solutions  match  well. 

As  can  be  seen  from  Figures  5.2,  5.3  and  more  from  Figure  5.4, 
the  oscillations  of  the  cross  section  as  a  function  of  ka  do  not  match 
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very  well.  Actually,  the  tendency  is  the  same  but  the  minima  and  maxima 
are  shifted  as  ka  increases  when  compared  with  the  GTD  solution  in  the 
backscatter  case,  Meixner's  solution  gives  a  much  better  result  for 
large  disks. 

For  a  disk  of  given  size,  ka,  the  agreement  of  the  E  plane  cross 
section  is  rather  good  even  for  large  values  of  ka  as  shown  in  Figure 
5.5,  5.6  and  5.7.  The  greater  difference  between  the  two  solutions  in 
this  E  plane  appears  to  be  the  value  of  the  minima. 

On  the  other  hand,  the  ^-components  do  not  match  very  well  for 
values  of  0  greater  than  20°  for  large  values  of  ka.  Figure  5^.8,  for 
ka=2,  showssgood  agreement  but.  Figure  5.10,  for  ka=10,  shows  H-plane 
values  obtained  from  Flammer's  solution  much  larger  than  those  obtained 
from  Meixner's  as  e  approaches  90°.  The  difference  almost  reaches  a 
factor  2.5  (i.e.,  ^s4  dB)  for  es=90° . 

We  have  seen  that  this  difference  appears  also  in  the  H-plane 
cross  section  for  ka=4  and  in  the  4>-part  of  Figure  5.12.  The  good 
agreement  for  ka=2  makes  less  likely  an  error  in  the  programming  of 
Flammers  solution,  but  nevertheless,  that  possibility  cannot  be  rejected. 
When  (f>  =0°  and  os=90°,  [E]^  is  proportional  to  ei-  An  error  due  to  a 
truncation  or  a  precision  problem  may  appear  in  the  computation  of 

3*  for  large  values  of  ka .  In  our  program  we  used  the  same  criterion 

to  truncate  the  infinite  summations  as  Hodge' s[9],  and,  therefore,  one 
would  not  expect  that  this  should  lead  to  a  difference. 

The  phases,  Figures  5.13  to  5.18,  seem  in  general  on  good  agree¬ 
ment  for  all  values  of  es,  except  for  an  evident  problem  in  Figures 
5.13,  5.14  and  5.18.  The  phase  of  the  3-component  of  the  field  obtained 
from  Meixner's  solution  is  consistently  180°  out  of  phase  with  that 
obtained  by  Flammer.  Thus  there  is  apparently  a  sign  problem  in  one  of 
the  components  of  Hodge's  since  he  does  not  obtain  the  same  phase  for 
the  backscatter  case,  es=0° ,  for  ^O0  and  <j>s=90°  as  can  be  seen  by 
comparing  Figures  5.13  and  5.15,  for  ka=2,  and  Figures  5.14  and  5.17 
for  ka=10.  We  verify  that  the  scattered  electric  field  must  have  the 
same  phase  for  <f>s=0°  and  4>s=90°  when  es=0°: 


.=0°  E5  =  Es  e  =  Es  e 


h  =90°  E5  =  E2  e  =  Es  e 

ps  o  e  o  y 


(129) 


Flammer's  solution,  as  calculated,  satisfies  this  condition. 

From  the  data,  we  see  that  we  have  an  overall  correspondence  be¬ 
tween  Flammer's  and  Meixner's  solutions  for  the  bistatic  normal 
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incidence,  parallel  polarization  case.  However,  except  in  the  case  of 
small  electric  circumference,  ka,  the  curves  are  never  quite  equal. 

The  weakest  point  of  the  solution  is  apparently  the  behavior  of  the  $- 
component  of  the  scattered  far-field  in  the  H  plane  in  the  neighborhood 
of  the  plane  of  the  disk. 


D  Conclusion 


In  this  part,  we  have  rederived  Hammer's  solution  for  the  normal 
incidence,  perpendicular  polarization,  bistatic  case.  We  have  shown 
that,  with  an  appropriate  choice  of  the  two  expansions  of  the  incident 
plane  wave,  the  condition  on  the  ^-component  of  the  total  electric  field 
at  the  edge  is  equivalent  to  the  whole  edge  condition.  This  remarkable 
property  will  not  appear  in  Meixner's  solution  as  we  will  see  in  the 
next  section.  It  is  a  characteristic  of  Hammer's  approach.  No  answer 
'oS  been  given  to  the  question  of  the  validity  of  the  general  solution 
for  arbitrary  incidence  and  polarization.  The  computational  checks 
explained  in  Section  B  should  allow  a  better  understanding  of  Hammer's 
approach.  The  numerical  test  of  Hammer's  electric  scattered  far- 
field  has  led  to  mixed  conclusions.  Despite  an  overall  agreement  with 
Meixner's  solution,  Flammer's  fields,  show  some  important  deviations. 
These  might  be  due  to  computer  problems. 

We  will  now  consider  Meixner's  considerably  different  approach 
to  the  problem  of  the  scattering  by  a  disk. 
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Figure  5.2.  Backscatter  normalized  cross  section  as  a 
function  of  ka. 
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Figure  5.4.  Normalized  cross  section  as  function  of  ka 
for  *s=90°;  es=45°. 
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Normalized  cross  section  as  function  of  es 
for  ka=4;  <j>s=90°  (E-plane). 
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igure  5.13.  Phase  of  Ej*  as  a  function  of  es  for  ka= 
4>s= 90°. 
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CHAPTER  VI 
MEIXNER'S  SOLUTION 


The  validity  of  Meixner's  solution  has  been  established.  Theo¬ 
retical  and  numerical  checks  have  been  conducted  in  recent  years,  and 
they  have  shown  very  good  agreement  between  computations  and  experi¬ 
mental  data.  In  this  Chapter,  we  will  not  repeat  Meixner's  solution, 
since  the  derivation  is  available  in  the  1 i terati^e.  The  reference 
work  used  for  Meixner's  solution  is  Hodge's  version[8].  Here  we  will 
derive  the  expansions  of  the  fields  in  terms  of  spheroidal  vector  wave 
functions  and  show  a  proof  of  the  solution  using  only  vector  wave  func¬ 
tions  without  refering  to  vector  potentials. 


A.  The  Vector  Potentials 


Meixner  used  Hertz  vector  potentials  to  solve  the  scattering  pro¬ 
blem.  In  this  section  we  will  give  the  expressions  of  the  components 
of  those  potentials  in  terms  of  the  scalar  wave  functions.  Equation  (7) 
This  formulation  is  much  more  convenient  than  Meixner's  for  the  compu¬ 
tation  of  the  fields.  Meixner's  solution  computes  the  potential  of  the 
scattered  field  for  arbitrary  incidence  and  polarization.  We  will  use 
here  the  notation  defined  in  Chapter  V-B ,  Equation  (114),  for  the  in¬ 
cident  field. 


E1  =  Eq (cost*  ElH  + 


rl  1  \ 

sinct  E  ) 


(130) 


The  definitions  of  E1  ,  Equation  (114),  in  this  work  and  in 
Hodge's  work[8]  have  opposite  signs.  Our  definition  was  chosen  to 
avoid  a  supplementary  minus  sign  in  the  expansion  of  E1"  in  vector 
wave  functions  in  Flammer's  solution.  This  will  lead  us  to  introduce 
a  minus  sign  in  front  of  the  scattered  vector  potentials  corresponding 
to  parallel  polarization.  As  in  Flammer’s  solution,  we  can  obtain  the 

-SI 


scattered  field  in  any  arbitrary  case  if  we  know  E  $nd  E^-1,  scattered 
fields  corresponding  to  the  incident  fields  E1"  and  E1j-,  respectively. 
We  will  define  the  vector  potentials  in  order  to  obtain  ESi.  directly. 


Meixner's  solution  involves  three  different  potentials,  n 
the  potential  of  ^he  incident  plane  wave,  tt  is  the  "reflected" 


-i 


is 


72 


i  s  ] 

potential  defined  so  that  n  +n  has  its  components  tangential  to  the 
disk  equal  to  zero.  Because  of  the  shape  of  ni,  the  components  of  the 
electric  field  corresponding  tangential  to  the  disk  to  iP+iT5  are  then 

equal  to  zero.  tts2  is  a  second  scattered  potential  defined  so  that 
7r=Tr ^  +ttS  1  +n s 2  satisfies  the  edge  condition  derived  by  Meixner: 


ft  <*p>  ■ 0  ■  s  <"»> at  <’3,i 

, 

ttx  and  7Ty  finite  at  the  edge. 

where  -rrp  =  irx  cos<i>  +  Try  sin<j>  .  j, 

•iP2  also  leads  to  an  electric  field  whose  components  tangential 
to  the  disk  are  equal  to  zero.  We  list  here  the  vector  potentials  for 

— c  II  j 

El  as  defined  above,  using  Hodge's  notation 


(132) 

(133) 

(134) 

(135) 

(136) 

(137) 


i 


73 


_  (  cn 

$2  o  v  .-m 

"•i  =  ”7  ).  1 

iy  k  m=0 


[Xm+l+(1‘5m,l)Xm-l]  cos  « 

.11  II  . 

1 _  !  Uaai  +  (1+*  i)^1  il  sin 

——  lum+l  m,l'  m-lj 


C6S-0 


(138) 


The  different  functions  used  in  those  potentials  are: 
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(144) 


As  previously  mentioned,  the  sign  of  n^,  Equation  (133),  is  the 
opposite  of  Hodge's  result  in  order  to  match  the  different  definitions 
of  EP" .  However,  the  signs  of  the  coefficients  of  and  have  not 
been  modified  in  Equations  (137)  and  (138)  in  order  to  keep  the  parallel¬ 
ism  between  the  parallel  and  perpendicular  cases.  The  sign  difference 
will  appear  inside  U*'  and  x"  instead.  We,  however,  have  the  same  Ll" 

i  mm  m 

and  X^  as  in  Hodge's  results.  This  is  due  to  a  sign  error  in  the  deri¬ 
vation  of  u"  in  Hodge's  solution  and  to  the  fact  that  the  x"  are  all 

m  3  ii  i  m  — q? 

zero.  In  order  to  simplify  U"  and  XA,  their  coefficients  in  tt  differ 

3  m  m3 

by  the  factors  (co/k  cos0Q)  and  (eQ/k  )  respectively  from  Hodge's  ex¬ 
pression.  We  now  can  express  these  potentials  in  terms  of  the  scalar 
spheroidal  wave  functions. 

From  Equations  (7),  (42)  and  (139)  we  can  express  Vm  cos  m4>  as  a 

summation  of  functions: 
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In  the  same  way,  from  Equations 


(7)  and  (140),  we  have  for 
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By  using  these  expansions  in  the  expressions  of  the  scattered 
potentials.  Equations  (133)  to  (138),  we  obtain  the  following  scattered 
vector  potentials. 


1 .  Perpendicular  polarization 


Si  _ 


i  ct> 

-T  l  i~m 

k  m=0 


«>,  -n 


^Xrn+r(1'6ml^Xrri-l^ 


,  j  S  (-ic  ,0)  mn 
n=m  Nmn  mn 


Rin)(-1c'i0>  (4) 

R(4>  (-1C.10)  J°"n 


(147) 


Si 


f.  » 


<2  Jo  n=m  Y""(  o’  R<4>HciO) 


(4) 

mn 


*7!  rm 

k 


m=0 

.  n 


cCi+(1-‘.i)xi-i] 


*  J  i— S  (-ic,0)  4V -  *e(4) 

n=m  NKn  ™n  R(4)(-ic,i0)  e,nn 

mn 


(148) 


2.  Parallel  polarization 
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Those  formulas  will  allow  us  to  compute  the  fields  simply  by 
taking  their  curls.  This  formulation  is  very  convenient  for  this  pur¬ 
pose  since  to  each  scalar  wave  function  in  a  given  component,  x  or  y, 
of  the  potential  corresponds  a  vector  wave  function  of  vector  ev  or  e  , 
respectively,  when  the  curl  of  the  vector  potential  is  taken.  The  ^ 
expressions  for  the  electric  and  magnetic  field  will  then  be  obtained 
in  a  very  straight  forward  manner. 


B.  Meixner’s  Fields 

_ ||  - 

In  this  section,  we  will  derive  the  expressions  of  II  and  E^A 
as  summations  of  vector  wave  function.  The  fields  are  computed  from 
the  Hertz  vector  potential  using  the  following  general  relations: 


H  =  i  ll)  V*TT 
E  =  —  VxVXTT 


(151). 


Using  Maxwell’s  equations,  we  also  have  the  relation: 


H  =  ft-  vxE  (152) 

o 

where  ZQ  is  the  free  space  wave  impedance. 

Equation  (152)  will  allow  us  to  compute  the  magnetic  fields  from  the 
trio  fields.  The  scattered  magnetic  fields  will  not  be  listed 
<■<  a  j-.<*  of  the  simplicity  of  the  transformation.  In  order  to  compute 
•'«  ‘  i>Js ,  we  will  use  the  following  identity: 
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x  or  y  . 


(153) 


The  final  form  of  the  expansion  of  the  fields  will  be  obtained  by  re¬ 
placing  the  vector  wave  functions  Ne*  and  Ny  by  their  expressions 
+  __  omn  omn 

in  terms  of  Np  ,,  and  Ne",  .  This  way  the  electric  field  will  be 
gm+ l n  Qm- i n 

expanded  in  terms  of  functions  whose  :-dependence  is  indicated  by  their 
m  index.  These  transformations  are  given  by  Flammer[10]: 
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We  will  now  list  E51  and  p\  obtained  by  applying  Equations  (151) 
and  (153)  to  the  corresponding  vector  potential,  respectively.  Equations 
(147)  and  (148)  and  Equations  (149)  and  (150). 


1  .  Perpendicular  polarization 
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We  notice  that  the  first  summation  is  a  "reflected" 
in  Chapter  V. 


field  as  defined 


Using  Equation  (154)  in  the  above  formula  and  rearranging  in 
terms  of  Xj,  we  obtain: 
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with  XA  =  0  for  m  £  0. 
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2.  Parallel  polarization 
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In  the  same  way,  using  Equation  (154)  in  the  above  formula  and 
rearranging  in  terms  of  U^,  we  obtain: 


csii  = 


kcosQ 


o  Jo  X  ci 


mn 


*  %-ln^ 


t  ,_L_  J 
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w*'  .h  U  =0  for  m  <  0. 
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3.  Conclusion 


For  an  incident,  field  of  arbitrary  incidence  and  polarization, 
the  total  scattered  E-field  can  be  expressed  in  terms  of  E^1  and  E3" 
as  shown  in  Equation  (115): 


Es  =  EQ(coSa  15s"  +  sin™  E^1) 

=  E  (cosa  H5"  +  Sina  H5") 

11  « 
where,  from  Equation  (152),  Hr1  =  (i/kZ  )  VxtL. 


(159) 


In  this  section  we  have  found  the  expression  for  the  scattered 
electromagnetic  field  in  terms  of  the  spheroidal  vector  wave  functions. 
This  allows  us  to  deal  with  the  fields  directly  instead  of  the  vector 
potentials,  and  it  therefore  offers  the  same  possibilities  as  Flammer's 
solution  for  studying  the  near-field  problem  in  a  much  simpler  manner. 

In  the  normal  incidence  case,  we  now  have  two  different  formulas, 
Flammer's  and  Meixner's  formal  solutions,  giving  the,  bistatic  scattered 
field.  A  term  by  term  comparison  of  the  two  solutions  should  be  carried 
out  to  see  how  it  is  possible  to  match  the  results.  A  first  step  how¬ 
ever  would  be  to  compare  the  scattered  electric  fields  in  the  far-field 
region.  This  would  be  much  easier  since  the  remaining  components  of 
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E5  are,  in  each  solution,  equal  to  single  summations  of  scalar  spher¬ 
oidal  functions  over  n.  The  task  of  reordering  those  far-field  com¬ 
ponents  as  summations  of  the  Legendre  associated  functions  has  not  yet 
been  completed.  The  Legendre  associated  functions  are  the  natural 
functions  since  they  are  orthogonal  in  the  far-fiela  region  and  since 
the  spheroidal  coordinate  system  approaches  the  spherical  coordinate 
system  in  this  limit.  This  orthogonality  will  lead  to  an  infinite 
number  of  relations  between  the  expansion  coefficients  dPn(-ic)  and 
values  of  the  radial  functions  at  c=0.  The  validity  of  these  relations 
could  easily  be  checked  on  the  computer  since  the  different  functions 
are  readily  available  in  Hodge's  program[9]  and  in  the  one  used  in 
Chapter  V  of  this  work.  A  study  of  those  relations  as  function  of  the 
variable  c=(kd/2)  might  explain  the  discrepancy  obtained  in  Figure  5.1 
for  the  backscattering  cross-section.  It  would  determine  whether  the 
differences  are  due  to  some  inaccuracy  in  the  computations  of  the  coef¬ 
ficients  or  to  improper  truncations  of  the  infinite  summation^  in  one 
of  the  solutions.  The  far-field  of  Meixner's  solution  can  easily  be 
computed  by  inserting  the  large  argument  approximations  of  the  spher¬ 
oidal  wave  functions.  Equation  (34),  in  the  formulas  for  these  fields. 
This  showed  a  sign  error  in  the  9-component  of  the  scattered  field  in 
the  perpendicular  polarization  case  as  computed  by  Hodge[9].  A  sign 
error  also  appears  in  the  e-component  of  the  scattered  field  in  the 
parallel  polarization  case  in  Hodge's  work. 

In  the  next  section,  we  will  derive  Meixner's  solution  with 
another  method  using  vector  wave  functions  exclusively.  The  edge 
condition  will  involve  all  the  components  of  the  electric  field  unlike 
Flammer's  solution  where  the  ^-component  only  was  needed. 


C.  Another  Proof  of  Meixner's  Solution 


As  can  be  seen  from  the  previous  section,  Meixner  develops  a 
solution  to  the  scattering  problem  by  using  only  one  expansion  for 
the  incident  plane  wave  in  each  polarization  case  -  parallel  or  per¬ 
pendicular.  His  method  is  equivalent  to  adding  to  the  incident  and 
reflected  field  a  summation  of  vector  wave  functions  whose  tangential 
components  are  zero  on  the  surface  of  the  disk  but  which  has  a  sin¬ 
gularity  at  the  edge.  The  variable  coefficients  of  that  summation 
are  then  adjusted  to  cancel  the  singularities  in  the  ^-component  of 
the  total  electric  and  magnetic  fields,  and  to  insure  a  proper  be¬ 
havior  of  their  other  components  near  the  edge.  This  is  exactly  the 
method  that  we  are  going  to  use  to  prove  Meixner's  solution. 

In  this  section  we  will  first  derive  some  useful  vector  wave 
functions  expansions  which  are  equal  to  zero  on  the  surface  of  the 
disk.  We  will  summarize  the  power  series  approximations  of  their 
components  in  terms  of  n  and  S  in  the  vicinity  of  the  edge.  We  will 
then  be  able  to  prove  Meixner's  solution,  using  the  above  expansions 
to  match  the  edge  condition. 
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1 .  Vector  wave  expansion 


This  section  is  mainly  based  on  the  properties  of  the  Nx  ex¬ 
pansion  of  a  plane  wave.  Equation  (47),  for  an  angle  of  inciS?Rce 
o  equal  to  90°. 

a.  Derivation  of  the  vectorial  relations 

— xM  1 

We  will  first  establish  some  dependence  relations  for  the  N  '  ' 

emn 

vector  wave  function  and  use  the  ^-orthogonal i ty  to  obtain  relations 
for  each  value  of  m.  We  will  then  show  that  those  relations  apply 
equally  to  the  odd  vector  wave  functions.  We  can  then  find  vector 
wave  expansions  whose  tangential  components  are  equal  to  zero  on  the 
surface  of  the  disk  and  that  behave  like  outward  travelling  waves 
in  the  far-field  region  by  using  a  term  by  term  matching  technique 
identical  to  that  used  for  the  calculation  of  the  reflected  field 
in  Flammer's  solution. 

For  o  =  90°,  Equation  (47)  multiplied  by  cose  can  be  written 
as  0 

r  V  ,  (-/?)  Nx^)(  %r  ,*)  =  5  f  1601 

p»-|  7  y 

»i  n  r»  r  "nn 


From  Equation  M?)  and  Equation  (14),  we  have,  for  n-m  odd, 
*mn("/2)=0  and,  therefore,  in  the  previous  summation  only  the  terms 
for  n-m  even  will  remain. 
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l  '  'mn(V2)  Nx(1)(n,C,F) 
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(161) 


We  will  introduce  here  the  functions  N. 
a  convenient  ^-dependence,  Equation  (154). 
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mn 


and  N 
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which  have 


Equation  (161)  can  be  rewritten  as  follows 
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(162) 
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l. 


.  *v: 
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where  corresponding  components  of  and  Nq^  have  the  same 

dependence. 


mn 


mn 


A  difficulty  encountered  with  the  use  of  this  notation  is  the 
difference  in  the  parity  of  n-m  between  Equation  (161)  and  (16?). 

Equation  (161)  is  a  summation  of  N*'  '  for  n-m  even,  while  Equation 

_ f.  ( i  \  i  \ 

(162)  is  a  summation  of  Nfi'  and  N~'  for  n-m  odd.  In  order  to 

avoid  a  possible  confusion^the  inde$np  will  be  used  instead  of  n 

in  N^^and  N"^  during  the  rest  of  this  section.  With  that  notation, 
8™  emn 

the  <f>- orthogonal ity  applied  to  Equation  (162)  leads  for  every  value 
of  m  to  the  relation: 


i  (1+6,  )( 1-6  )y  ,  (rr/2)N 
=m-l  lm  om'  m-lpl  '  '  i 


p=m-l 


;0) 

"mp 


p=m+1 


Ym+lp 


(V2)N  (1)  =  0 
mp 


063) 


Equation  (163)  means  that  the  coefficient  of  cos  in  the  n- 
and  ^-components  of  the  above  summation  and  the  coefficient  of  sin 
m$,  except  for  m=0,  in  its  ^-component  are  equal  to  zero. 

We  can  then  substitute,  for  m  greater  than  zero,  sin  me  for  cos 
m 4>  and  cos  for  -sin  m<$>  without  changing  the  value  of  the  summation. 
This  transformation  leads  to  the  following  equation  for  the  odd  func¬ 
tions: 


0) 

mp 


+ 


p4»i  Tm+1p 


(»/2)N  (1)  -  ff 
mp 


(164) 


for  every  value  of  m£l  and  (n>5»4>)*  Note  that  Equations  (163)  and 

(164)  are  also  valid  for  This  is  proved  by  simply  taking 

gmp 

the  curl  of  these  equations.  In  order  to  match  the  behavior  of  out- 
wave  travelling  waves  at  infinity  and  still  keep  the  behavior  of  the 
tangential  components  at  f=0,  we  will  make  a  term  by  term  transformation. 
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r r om  Equation  (164)  and  Equation  (30)  to  1331,  we  ran  see  that,  on 

the  surface  of  the  disk,  the  n-  and  ^-components  of  fT+,’j  or  fT^’l 

§mP, 


are  proportional  to  R 
will  therefore  make 


Ri-1 (‘1C» io) »  or  Rm+  j  ( ' 1  c  >  i 0 )  respectively^  "We 
tne  following  substitution  in  Equation  (164): 
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We  will  denote  by  F  (  n,  6,  0)  the  value  of  the  transformed  sum¬ 
mation.  From  Equations  rl64)  and  (165),  we  have  for  every  value  of 
n: 
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F  satisfies  for  j  n  I  -  0: 

gm 


C^e  ]  (n,O,0 )  =  0 
0 

where  ct  =  n  ,0  .  ( 167) 

[F  ]  (n,O,0)  =  0  for  mi  1. 
gm  a 


In  contrast  to  Equation  (161)  and  (163)  these  F  functions 

gm 

are  not  zero  everywhere.  Note  that  F  is  the  reflected  field  cor- 

gm 

.  .-sponding  to  the  zero  incident  field  of  Equation  (164).  Since  the 
•fleeted  field  is  not  zero,  we  immediately  realize  that  the  reflected 
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NeM  is  not  equal  to  the  scattered  field  in  that  case.  We  will  need 
the  expansion  of  the  components  of  F  in  terms  of  o  and  £.  near  the 

§m 

edge.  In  the  following  subsection,  we  list  the  expansion  of  the  nec¬ 
essary  vector  wave  functions. 

b.  fT^^andiT^  near  the  edge 
8mn  8m- 1  n 


We  do  not  need  very  complete  expansions  of  the  vector  wave  func¬ 
tions  in  the  neighborhood  of  the  edge  since  we  are  only  looking  for  the 
terms  that  do  not  have  the  behavior  required  by  the  edge  condition 
for  each  component.  The  other  terms  will  not  lead  to  any  problem 
in  satisfying  the  edge  condition.  We  will  first  consider  the  case 

of  ,  component  by  component.  N^^has  no  singularity.  Since 

am+1n 

Equation  (161)  involves  the  functions  'for  n-m  even  only,  we  will 

_+  ( i )  _ t  ,•  \  mn 

need  N„  '  and  N  '  1  for  n-m  even  only  for  the  computation  of 
8m+ln  8m  -In 

F8m* 

Fo£  the  n-component,  we  only  need  the  terms  of  singularity  greater 
than  s-2,  where  s  is  the  distance  from  the  edge  as  defined  in  Equation 
(89). 


n-m  evenjnl  <<  1  and  cs  «  1 
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(168) 


gm+ln  n  ckd 


For  the  S -component  of  ,  we  only  need  the  terms  of  singu¬ 

larity  greater  than  s~ 2  om+^n 


ity  greater 
n-n  evenjnl  <<  1  and  cC  <<  1 


rn+(  i )  i  .  2i  I  cos,  J 

V’"1'  ‘  ckd2  J 


Smn(-ic’°)  C2-n2 


'+  0(  s'*5) 
(169) 


For  the  ^-component,  we  need  the  terms  up  to  the  order  of  s*5 
to  test  the  edge  condition  on  the  ^-component. 
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F  or  '-m  even ,  r; 
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Wp  will  need  the  expansion  of  ],  for  n-m  odd 
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(171) 


It  is  well  behaved  at  the  edge  since  it  goes  to  zero  as  s  does.  The 

case  of  can  easily  be.  derived  from  that  of  .  N™^ 

gm-ln  gm+ln  gm-ln 

has  no  singularity.  The  n-  and  S-components  of  TT^  can  be  obtained 

gm-ln 

from  Equations  (168)  and  (169)  by  replacing  m+1  by  m-1  in  the  -de¬ 
pendence.  The  expansions  of  the  '('-component  of  N~  are  obtained 

gm-ln 

from  Equation  M70)  and  (171)  by  replacing  m  by  -m  except  in  the  in¬ 
dexes,  s i n ( m+ 1 ) ❖  by  - s i n ( m- 1 ) and  -cosfm+l)^  by  cos(m-l)<>. 


In  order  to  compute  the  complete  electromagnetic  field,  we  will 

also  need  the  functions  M+fi)  and  M"^  as  defined  below: 

gm+ln  gm-ln 


As  seen  from  their  expressions  in  Flammer's  bookLlO],  the  n- 
and  ^-cgmponents  of  these  functions  have  a  singularity  of  order  at 
most  s” 2  and  therefore  satisfy  the  edge  condition.  The  expansion 

of  <t>-component  of  M4^  is  listed  below  for  n-m  even, 
gm+ln 
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We  will  also  need  W4^  L  for  n-m  odd: 

gm+ln  41 


n-m  odd 


^!!nVe  Smn ^ " 1  c * 0 ^  W'ic’io)  js?>+1)' 


r=  0(s)  (174) 


The  corresponding  expansions  for  are  obtained  by  replacing 

gm-ln 

m+1  by  m-1  in  the  trigonometric  functions.  We  note  that  the  components 
of  the  functions  of  index  4  listed  in  Equations  (168)  to  (170)  and 
(173)  have  singularities  that  do  not  satisfy  the  edge  condition,  but 
that  all  the  components  of  the  vector  wave  functions  of  index  1  do 
satisfy  it. 


F4  is  expanded  in  terms  of  N_^ instead  of  N4^  or  N~^  . 

8m  gmp  gm+ln  gm-ln 

From  N4^  ,  n-m  even  or  odd,  we  can  derive  the  expression  of  N4^ 
Bm+ln  gmp 

for  p-m  odd  or  even,  respectively,  by  replacing  m  by  in- 1  even  in  the 
indexes  and  n  by  p.  From  TT'  n-m  even  or  odd,  we  can  derive  the 

-m  8m_1 

expansion  of  N  '  'for  p-m  odd  or  even,  respectively,  by  replacing 
gmp 

m  by  m+1  and  n  by  p. 
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c.  F  at  the  edge 

8tn 

Here  we  will  compute  the  coefficients  of  the  improper  singularities 
of  the  components  of  F  and  1/kv  x  F  .  The  coefficient  of 

rK;  Bm 

— 'T'cTo  in  [F  ]  is,  from  Equation  (168)  and  (166): 
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the  coefficient  of  — o‘^5 -t/o'  ’n  *-  E0  -i-  is  halr  of  equation  1/5;. 
(n^)b/^  §m  ’ 

Those  two  components  will,  therefore,  enforce  the  same  condition. 

The  coefficient  of  in  [E  ]  is,  from  Equation  (166)  and  (170): 
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It  can  be  rewritten  as  follows: 
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This  coefficient  has  a  different  ^dependence  than  the  other  two. 

The  coefficient  of  — =■  ~  in  Cr  E  1,  is  obtained  from  Equation 

n  +;.  K  §m  v 

(166),  (172)  and  (173).  It  is  equal  to  Equation  (175)  multiplied 
ky  2  ‘  Equation  (175)  and  (177)  differ  only  by  their  ^-dependence 
and  a  factor  m  in  Equation  (177). 

Equation  (175J_  and  (177)  give_the  coefficients  of  the  terms  in 
the  components  of  F  and  1/k  v  x  F  that  do  not  behave  as  required 
om  om 

by  the  edge  condition  in  the  vicinity  of  the  edge.  Note  that  for 
m=0,  the  coefficient  in  Equation  (177)  is  always  equal  to  zero  while 
only  the  odd  case  of  Equation  (175)  disappears.  Thus  Fn  and  1/k 
V.x  Foo  behave  as  required  by  the  edge  condition  and,  therefore,  the 
missing  relation  in  Equation  (167)  will  not  be  needed  since  the  func¬ 
tions  TT  and  M  are  well  behaved  at  the  edge, 
op  op 

We  now  have  the  necessary  tools  to  solve  the  scattering  problem. 


2.  Solution 


The  solution  will  be  derived  for  the  general  case  of  arbitrary 
incidence  and  polarization.  We  will  use  the  same  notation  as  in  part 
V  for  the  electric  and  magnetic  field. 

E1  =  Eo(COSa  E1"  +  Sina  E1a)  (178) 

We  expand  E1  and  f1  from  Equations  (47)  and  (48) 
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From  Equations  (30)  to  (33),  we  see  that,  for  n-m  odd,  the  in¬ 
cident  vector  wave  functions  have  null  taggential  components  on  the 
disk.  Therefore  the  reflected  fields  are: 
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i  V 


R^ih'ic.lo) 
mn  '  '  (jj-y(4) 


k  m=0  n=m  '*<"°)  ldJ>(-1c,1o) 


“mn 


(180) 
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Once  again  the  "reflected"  fields  as  defined  here  are  not  equ ; . 
to  the  scattered  fields  since  they  do  not  satisfy  the  edge  condition. 
We  can  now  add  to  the  reflected  field  any  summation  of  functions  F 
and  still  have  a  reflected  field.  The  coefficients  of  that  sum-  8m 
mation  will  be  adjusted  so  that  the  fields  satisfy  the  edge  condition. 


_ x  f  i  1 

From  Equation  (154),  we  see  that  N*'  '  involves  only  even  func¬ 
tions  N+^  and  N~^  while  1  ^  invokes  only  odd  functions 
em+ 1 n  em- 1 n  emn 


N*(l)  and  . 

m+ln  °m+ln 


Therefore  to  match  the  edge  condition,  we  need 

satisfies 


only  add  a  summation  of  F  to  Fr  and  F  to  Fr  1  sinrp  r 

em  om  c  »  since  r 

gm 

~ne  scattered 


the  boundary  condition  on  the  surface  of  the  disk. 


;  e  ■  l  h  ' 


o  e: 


ES" 


-r" 


i 


k  cos o. 


V 

m=0 


(181) 


=  E' 


k  J.  X  F 

m=l  mv  o 


m 

where  Umv  and  Xmv  are  adjusted  to  satisfy  the  edge  condition. 

The  summation  over  m  in  FSl  begins  with  m=l  since  F  is  well 

°o 

behaved  at  the  edge  and  does  not  correct  the  behavior  of  Ert  We  will 
successively  consider  the  perpendicular  polarization  and  parallel 
polarization  cases. 


a.  Perpendicular  polarization 

Equations  M  7  9 )  and  (1801  can  be  rewritten  as  follows  from  Equ 
♦hvi  1541:  1 
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(182) 


We  must  adjust  the  coefficient  X  so  that  the  coefficients  of 
the  ill-behaved  terms  in  each  component  disappear  at  the  edge. 

Taking  advantage  of  the  ^-orthogonal ity,  we  can  write  a  specific 
condition  for  each  value  of  m: 

for  every  m  10,  atn  =£=  0: 


1 

^  p=m-l  0+6imHl-iom)'rm_1p(e0)  N0^ 

-r  LWV^'1 

p=m+l  r  mp 


iilL(-ic,io)  , 
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,  Rm-  o("ic’io)  _ (4) 

'  *  p=Ll  [(Wta)(M°">)Weo)  "“mp 


+  1  X  f  (1-  <5  )  must  satisfy  the  edge  condition. 
k  mv  om  '  om' 
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We  have  seen  that  the  functions  of  index  i  and  behave  as 

°op 

required  by  the  edge  condition.  The  absence  of  coefficient  X  is 
therefore  not  a  problem.  Tor  any  other  value  m,  m>l,  the  conations 
on  X.k  from  the  three  components  ’W"  of  tquation'7  183)  are  identical 
.me"  Me  coefficients  of  the  improper  singularities  in  Equation  (1681 
to  (I  ' * )  are  proportional  apart  from  'heir  ^-dependence.  From  Equation 
(173j,  we  see  that  this  also  applies  for  the  o-component  of  the  curl 
of  Equation  (183)  which  corresponds  to  the  m-part  of  the  expansion 
of  the  magnetic  field. 


Xmv  satisfies,  for  mil: 


1  Sm-lp(-tc-°) 

"Sp<-1c"°> 

oo 

*  pi,  WV  1 
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p=m_l  (1  +  '  lm^^'om^m-lp^o^  c  ^[47 
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+  xmvi  pj‘,  f1+go-gviP(^)  1 


p-  n+  i 


c  m+lp 


Rm+ 1 p ■ " 1 c ’ 1 0 ' 


Rm-lp(-ic’i0' 
=  0 


We  notice  that,  from  Equation  (42),  we  have,  m  ._0: 

n-OO-^ijY™  ■,  ( e  )  =  J1”  Sm  ,  (-ic.cose  )(l-s  ) 
v  om7  lm7  m-!pv  o  N  ,  m-ln  *  o  om 
r  m-lp  K 

'm+,p<  6o>  '  «£—  S»HpHc,C0S80) 


(185) 


Introducing  the  functions  W  and  as  in  Meixner's  solution, 
Equations  (1431  and  (144),  we  have  for  ”  : 


mv 


Wm-1  "  Wm+1 
^m-l  +  ^m+1 


for  m  >  1 


(186) 


Witn  those  coefficients  we  can  compute  the  scattered  field  of 
f  i"  jn'qvii  In  .jlir  polarization  case.  Wo  obtain  from  Equations  (180), 
f  181 ;  and  ( Mfi) : 
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(187) 

We  can  check  this  solution  by  comparing  it  with  Equation  (156). 

We  compare  the  coefficients  of  and  fT^  in  both  solutions. 
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The  scattered  field  obtained  by  this  method  is  identical  to  the 
one  obtained  by  Meixner,  as  expected. 


b.  Parallel  polarization 

Equations  (179)  and  (180)  can  be  rewritten  as  follows,  from  E- 
quation  (154): 
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We  must  adjust  the  coefficients  IJ  in  Equation  (1ST)  so  that 
toe  coefficients  of  the  improper  singuT^ri t ies  in  each  component  of 
the  fields  disappear  at  the  edge. 

Taking  advantage  of  the  '{‘-orthogonality,  we  can  write  a  specific 
condition  for  each  value  of  m: 


tor  m  t)  at 


mu 
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mp 


--  U  F 


k  cosen  mv  em 


si  satisfy  the  edge  condition. 


the  functions  of  index  1  behave  as  required  by  the  edge  condi¬ 
tion  and  will,  therefore,  not  contribute  to  the  computation  of  U  . 

For  m  greater  than  zero,  the  conditions  on  U  from  the  three  com¬ 
ponents  of  Equation  (190)  and  from  those  of  its  curl  are  equivalent 
since  the  coefficients  of  thq  improper  singularities  in  Equations  (168) 
to  (170)  are  proportional  apart  from  their  ^-dependence.  For  m  equal 
to  zero,  the  n-  and  C-components  of  Equation  (190)  and  the  ^-component 
of  its  curl  lead  to  equivalent  conditions  on  (I  .  The  ^-component 
of  Equation  (190),  however,  contains  no  improper  singularity  since 


the  transformation  of  the  coefficient  of  the  singularity  of[<(4)  1  , 

__(4)  m-l,n  9 

Equation  (170),  into  that  of  Cn^  3, leads  to  a  m  factor  in  the  latter. 

mp  9 

It  can  be  seen  in  [F  ]  ,  Equation  (177),  for  example.  For  m  equal 

emp 

to  zero,  tiie  ,|, -component  of  Equation  (190)  does  not  impose  any  con¬ 
dition  on  U  in  order  to  satisfy  the  edge  condition.  The  satisfaction 
of  the  threeother  conditions  will  then  be  sufficient  at  the  edge 
to  obtain  the  scattered  field.  For  any  value  of  m,  the  condition 

on  !l  i 
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Introducing  the  functions  W  and  <|<  as  in  Meixner's  solution. 
Equation  (143)  and  (144),  and  using  the  properties  of  Ymn(eo),  Equa¬ 


tion  (185),  we  have  for  Umv: 

W1 
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ov  4,1 

W  ,  +  W  , 

U_..  =  JUL.- m-tl.  for  m  >  1 
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(192) 


With  those  coefficients  we  can  compute  the  scattered  field  of 
the  parallel  polarization  case.  We  obtain  from  Equations  (180),  (181) 
and  (192): 
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We  can  check  this  solution  by  comparing  it  with  Equation  (158). 
We  compare  the  coefficients  of  tC(4)  and  r(4)  in  both  solutions. 
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(194) 

We  see  that  the  coefficients  match  exactly.  The  scattered  field 
obtained  by  this  method  is  exactly  that  of  Meixner's  solution. 

We  have  given  a  proof  of  Meixner's  solution  by  dealing  with  electro¬ 
magnetic  fields  only.  The  scattered  field  has  been  determined  so 
that  it  satisfies  the  boundary  conditions  on  the  surface  of  the  disk 
and  the  edge  condition.  This  proof  gives  more  insight  in  the  behavior 
of  thp  fields  in  the  vicinity  of  the  edge  than  Meixner's.  It,  however, 
requires  the  calculation  of  the  expansions  of  the  spheroidal  vector 
wave  functions  in  the  vicinity  of  the  edge,  and  its  edge  condition 
consists  of  six  conditions  whose  compatibility  must  be  checked.  It 
is  algebraical ly  more  complicated  than  Meixner's,  whose  edge  condition 
consists  in  one  equation  only  since  Meixner  mostly  deals  with  scalar 
functions.  In  the  next  section,  we  will  calculate  the  scattered  field 
for-  n  incident,  field  equal  to  a  single  vector  wave  function. 


Scattered  field  of 


jle  vector  wave  function 


From  the  solution  derived  in  the  previous  section,  we  can  easily 
calculate  the  scattered  field  for  an  incident  electric  field  equal 
to  a  single  wave  function.  The  purpose  of  this  section  is  to  show 
that  this  scattered  field  contains  all  the  vector  wave  functions  of 
index  4  which  have  the  same  index  m  in  their  t>- dependences  as  the 
incident  one.  Let  us,  for  example,  consider  an  incident  field  equal 

to .  Adjusting  Equations  (179),  (180)  and  (181)  to  our  case, 
K  em+ln 
we  have: 


r  r  = 


i  jC(4>  ,  Iki  f 

k  Rin^"ic» io)  em+ln  k  em+l 
ron 


(195) 


Following  the  same  derivation  as  in  the  parallel  polarization 
c  .tso ,  we  obtain,  from  Equation  (190),  the  following  condition  on  U  , 
fee  '  i  t  he  edge  emd  i  .  ion. 
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(197) 


From  Equation  (195)  and  the  definition  of  Fg  ,  Equation  (166),  we 

see  that  the  obtained  scattered  field  involv?s^all  the  functions 

<(*>  and  N  ^  for  the  given  value  of  m.  If  we  choose  an  in- 
gm+ln  gm+ln  _ 

cident  field  equal  to  r  N*'  ,  instead,  the  scattered  field  contains 

k  emn 

all  the  functions  and  fT^  for  m  given,  from  Equations  (154) 

em±  1  n  em±  1  n 

and  (195).  Whatever  the  vector  wave  function  of  vector  e  o re, 
or  derived  from  these,  Equation  (154),  the  corresponding  icatte^ed 
field  is  not  equal  to  the  reflected  field  and  contains  the  other  vec¬ 
tor  wave  functions  with  the  same  4>- dependence. 

D.  Conclusion 

In  this  part,  we  have  shown  that  the  fields  of  Meixner's  solu¬ 
tion  can  be  expressed  in  terms  of  spheroidal  vector  wave  functions. 

We  can  therefore  calculate  the  fields  everywhere  directly.  In  par¬ 
ticular,  this  will  allow  us  to  compute  the  near-fields  and  the  cur¬ 
rent  distribution  on  the  surface  of  the  disk.  Meixner  s  solution 
is  formally  as  convenient  as  Flammer's  in  this  region.  However,  for 
computational  purposes,  Meixner's  solution  sh°u1d^e  retained  since 
it  only  requires  the  values  of  the  angular  and  radial  spheroidal  func¬ 
tions  and  their  derivatives  for  n-m  even.  On  the  surface  of  the  disk, 
we  can  see  from  Equations  (23)  and  (24)  that  we  only  need  the  values 
of  the  radial  functions  at  5*0,  since  the  derivatives  can  be  expressed 
from  them.  We  have  also  shown  how  the  scattered  field  of  a  single 
incident  vector  wave  function  involves  all  the  vector  wave  functions 
of  the  same  ^-dependence.  This  is  a  consequence  of  the  fact  that 
the  vector  wave  functions  cannot  be  considered  to  be  modes  because 
of  the  dependence  relationships. 
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CHAPTER  VII 
CONCLUSION 


Considerable  knowledge  has  been  obtained  about  scattering  by 
a  metallic  disk  in  the  last  three  decades.  However,  most  of  this 
work  was  done  in  the  far-field  region. 

Meixner's  solution  can  be  directly  applied  in  this  case  because 
of  its  vector  potential  formulation.  The  purpose  of  this  work  was 
to  establish  a  valid  solution  everywhere  in  space. 

Flammer's  solution  was  considered  first  since  it  leads  directly 
to  formulations  of  the  scattered  electric  and  magnetic  fields.  The 
validity  of  this  solution,  however,  still  needed  to  be  established. 

We  rederived  it  for  the  normal  incidence,  perpendicular  polarization 
bistatic  case  and  compared  numerical  results  to  data  obtained  from 
Meixner's  solution.  The  general  agreement  of  both  solutions  showed 
that  Flammer's  solution  is  acceptable  in  this  case.  We,  however, 
encountered  problems  in  the  derivation  of  Flammer's  solution  in  the 
arbitrary  incidence  case  and  its  validity  is  still  questionable. 

Thus,  Flammer's  solution  did  not  provide  us  the  expressions  that  we 
were  looking  for. 

Meixner's  approach  to  the  problem  leads  to  a  solution  that  has 
been,  in  turn,  expanded  in  spheroidal  vector  wave  functions.  Another 
proof  has  been  given  directly  using  dependence  properties  of  the  sphe¬ 
roidal  vector  wave  functions.  This  formulation  of  the  fields  is  ade¬ 
quate  even  for  near  field  computations.  For  example,  the  computation 
of  the  surface  currents  on  the  disk  should  be  made  with  Meixner's 
solution. 

In  parallel  to  the  solution  of  the  scattering  by  a  disk,  we  have 
shown  some  dependence  properties  of  the  spheroidal  functions.  The 
derivatives  of  the  spheroidal  angular  functions  are  not  independent 
as  the  original  functions  are.  Some  dependence  relations  of  the  sphe¬ 
roidal  vector  wave  functions  of  index  1  have  been  derived,  and  it 
is  therefore  impossible  to  consider  a  single  vector  wave  function 
as  a  mode  in  these  cases. 

Some  questions  remain  after  this  work.  The  formal  agreement 
of  the  scattered  far-field  of  Flammer's  and  Meixner's  solution  in 
the  normal  incidence,  perpendicular  polarization  bistatic  case  should 


be  numerically  checked  as  explained  in  Section  VIB.  The  validity 
of  Flammer's  arbitrary  incidence  solution  should  be  tested  as  explained 
in  Section  VB.  If  the  validity  of  Flammer's  solution  can  be  established, 
it  should  be  interesting  to  compute  the  surface  current  distribution 
on  the  disk  from  Flammer's  and  Meixner's  solution.  This  could  be 
a  good  test  since  the  general  formal  solutions  are  too  complicated 
to  be  matched  term  by  term. 


99 


REFERENCES 


1.  G.  Mie,  Ann.  Physik,  25,  377  (1908). 

2.  Lord  Rayleigh,  Phil.  Mag.  44,  28  (1897). 

3.  J.  Meixner  and  W.  Andrejewski,  Ann.  Physik  7,  157  (1950). 

4.  J.  Meixner,  Ann.  Physik  (6),  6,  2  (1949). 

5.  C.  J.  Bouwkamp,  Philips  Research  Rep.  5,  401  (1950). 

6.  0.  B.  Hodge,  "Eigenvalues  and  Eigenfunctions  of  the  Spheroidal 
Wave  Equation,"  J.  Mathematical  Physics,  Vol.  11,  p.  2308  (1970). 

7.  C.  Flammer,  "The  Vector  Wave  Function  Solution  of  the  Diffraction 
of  Electromagnetic  Waves  by  Circular  Disks  and  Apertures,"  J. 

Appl .  Phys.,  Vol.  24,  pp.  1218-1231  (1953). 

8.  D.  B.  Hodge,  "Spectral  and  Transient  Response  of  a  Circular  Disk 
to  Plane  Electromagnetic  Waves,"  IEEE  Trans.  Vol.  AP-19,  p.  558 
(1971). 

9.  D.  B.  Hodge,  "The  Calculation  of  Far  Field  Scattering  by  a  Cir¬ 
cular  Metallic  Disk,"  Report  710816-2,  February  1979,  The  Ohio 
State  University  ElectroScience  Laboratory,  Department  of  Electrical 
Engineering;  prepared  under  Contract  N00014-78-C-0049,  Department 

of  the  Navy,  Office  of  Naval  Research,  Arlington,  Virginia. 

10.  C.  Flammer,  Spheroidal  Wave  Functions,  Stanford  University 

Press  (1957). - 

11.  J.  A.  Stratton,  Electromagnetic  Theory,  McGraw-Hill  Book  Company, 
Inc.,  New  York,  pp.  38-51  (1941). 

12.  C.  Flammer,  "The  Vector  Wave  Solution  of  the  Diffraction  of  E- 
lectromagnetic  Waves  by  Circular  Disks  and  Apertures,"  Report 
No.  27,  SRI  Project  No.  591,  September  1952,  Stanford  Research 
Institute;  prepared  under  Contract  No.  AF19(604)-266,  Air  Force. 


100 


APPENDIX  A 

THE  OBLATE  SPHEROIDAL  COORDINATE  SYSTEM 


Figure  A-l.  The  oblate  coordinate  system. 


The  oblate  spheroidal  cc  'dinate  system  has  the  components  (n, 
)  whose  characteristics  are: 

-  the  ^-constant  surfaces  are  ellipsoids  of  revolution  (around 
the  z  axis,  in  our  case).  The  intersection  of  those  surfaces  with 
a  plane  containing  the  z-axis  are  ellipses  whose  foci  are  the  edge 
points  of  the  circular  disk  in  this  plane. 
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-  the  'i-constant  surfaces  are  hyperboloids  of  re.olution  around 
the  z  axis.  In  the  same  way,  the  foci  of  the  intersections  with  a 
radial  plane  are  the  edge  points  of  the  circular  disk. 

-  :  is  the  usual  azimuthal  coordinate  as  in  the  cylindrical  case. 

;=0  represents  the  surface  of  the  disk  of  radius  d/2,  c-h  kd. 

varies  from  0  to  infinity  whereas  n  varies  from  -1  to  1. 

The  following  formulas  have  been  rederived  from  the  general  or¬ 
thogonal  coordinates  systems  transformations  as  explained,  for  example, 
by  St ra t ton [ill. 

We  first  establish  the  relations  between  the  usual  coordinates 
systems  and  the  spheroidal  coordinate  system..  We  then  consider  the 
limiting  cases  of  large  radial  argument,  c£,  of  the  surface  of  the 
disk,  £=0,  and  of  the  neighborhood  of  the  edge,  T=0  and  n=0. 

1.  Coordinate  Transformations 


a.  Components 
from  Cartesian: 

;  x  =  |  tr(i--r)(in2)f2  cos  $ 

t  y  =  f  [( 1- r>2) ( 1  +  C2)J%  sin* 
i  z  =  d  ,  f 

from  cyl indrical : 


(AT) 


(A2 ) 


from  spherical: 


r 


,  (V+z?  .  i 


j  l-Ai* 


cosm  =  n 


.  in- 


II-Cku-2)  1 

9  O  1 


,  2  2 
1-n  +•-. 


(A?) 


The  metric-coefficients  are  computed  from  the  following  formula: 


r 


h. 

1 


+ 


2  2 

+f— 5 ■  )  where  v^n,  £,  <P 
dv  J  \  3  v  / 


From  Equation  (Al),  we  obtain 


fn2t52  1 

.  d 

|V+«2  1 

1 

C\J 

cr 

i 

^  "  2 

i - 

PO 

i - 

[h«t>  =  |  f(l-n2)(l+C2)]^ 


(A4) 


(A5) 


b.  Unit  basis  vectors 


We  give  in  Table  A6  the  dot  products  between  the  basis  vectors 
of  the  spheroidal  coordinate  system  and  those  of  the  common  coordinate 
systems.  The  different  formulas  can  be  obtained  using  the  general 
relationship  in  a  basis  (e^): 

t  =  l  (a.e!.)e.  (A7) 

The  dot  products  are  computed  using  the  relation 


where  «  and  ?>  can  oe  x,  y,  z;  p,  r,  e,  n,  C,  - 
2 .  Large  c  Approximations 

a.  Components 

Since  £  is  always  positive  Equations  (Al)  to  (A3)  reduce  to: 


r  d 

.  2 

x  2 

«  J 

1-n  cos4> 

2 

1-n  s  i  n 

d 

p  =  ? 

cJ 

1-n2 

(A9) 

n  =  cosO 

sine  =  J  1-n^  for  cC  >;>  1 
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This  leads  us  back  to  the  spherical  coordinate  system.  The  two 
main  relationships  are: 

n  =  cose  and  £  =  (A10) 


b.  Basis  vectors 


For  C£  >>  1 

fe 


1* 


-nep  +  (1-n  )ez  +  -eQ 


(All) 


%  ->  er. 

On  the  Surface  of  the  Disk 


The  disk  is  represented  by  £  =  0.  We  have  the  relationships: 

p  =  r 


=  4il-n2 


e  >  - 

n 


-A-  0 
TnT  6p 


(A12) 


4.  Edge  Coordinate  System 


To  determine  how  a  component  of  the  field  behaves  in  the  vicinity 
of  the  edge  of  the  disk,  the  significant  quantity  of  interest  is  the 
distance  from  the  edge  to  the  observation  point.  Let  us  call  this 
quantity  s  in  the  coordinate  system  shown  in  Figure  A-2. 


z 
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We  will  use  the  following  notation: 


s  distance  from  the  edge  to  the  observation  point. 

t  angle  between  the  positive  side  (n>0)  of  the  disk  to  the 
observation  point. 

We  restrict  our  domain  to  small  distances  from  the  edge:  s  <<  d.  We 
notice  that,  in  this  system,  with  s  <<d: 

t  =  0  is  the  disk  itself 


t  =  ir  is  the  plane  of  the  disk  outside  the  disk. 

We  will  derive  the  relations  between  the  spheroidal  coordinate 
of  small  argument  and  this  edge  system. 

In  the  following,  we  have  C  <<  1 ,  j n |  <<  1 . 

By  definition,  we  have 


2  (  d  Y"  2 

s  =  |  p-  2  )  +z 


(A13) 


In  spheroidal  coordinates,  this  equation  can  be  written 
s2  =  j-  |[(l-n2)l5(HC2)Js  -  if  +  n2C2j 
For  n  and  (,  small,  we  may  write: 

fn-n2)(l  +  C2)]a  =  1  -  f-  +  f  +  0(n4,C4,n2,C2) 


which  leads  to 

2  _  d2  jl  /  ,2  „2\2  2  A  _  d2  ,  2if2v2 


(A14) 


(A15) 


s  = 


i  (£  "n  )  +n  **  f  =  ^6  (n  ) 


Final ly 


s  =  ^  (  n  +  f.  ) 


( A 1 6 ) 
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Gradient 

v  <t>  = 

Divergence 
v-F  = 

Curl 

vxF  = 


APPENDIX  B 

DIFFERENTIAL  OPERATORS  IN  THE  OBLATE  SPHEROIDAL 
COORDINATE  SYSTEM 
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APPENDIX  C 

ERRATA  IN  FLAMMER'S  WORKS 


The  errata  contained  in  this  appendix  have  been  collected  by 
Dr.  Hodge,  Dr.  Garbacz  of  the  ElectroScience  Laboratory  and  myself. 
The  following  corrections  should  be  made. 

A.  Errata  in  Flammer's  Book  [10] 


Page  42  -  Equation  (4-6-14)  should  be  replaced  by  the  following,  for 
n-m  odd: 


•n-m-1  JB  .  m+1  .mn  ,  .  , 

(1),  .  1  2  «"'•  c  di  (~ic) 

mn  v  ’  '  <»,  .  .  * 


(2m+3)  l  '  ifftU 


(Cl) 


r=l 


and  Equation  (4-6. 15b)  by,  for  n-m  odd:  » 

1n-"+,(2m+3)  1‘  d"n(-ic) 


R^’Mc.lo)  .  1 


r=l 


cR^^-ic.io)  2m  m!  cm+2  d™n(-ic) 


(C2) 


Page  43  -  ji  should  be  replaced  by 

n-m- l-  j |  -n  RU)  for  n_m  0dd,  Equation  (4.6.16b) 

Page  47  -  The  ^-dependence  term  cos  m(  4>-  <t>' )  should  be  introduced  in 
the  expansions  of  the  Green's  functions  in  Equation  (5-2-11)  and 
(5-2-12). 
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B.  Errata  in  Flammer's  Paper[7] 

Part  I 

Page  1219  -  Read  ir+m_n  instead  of  (-l)r+m_n  in  the  definition  of 
Rmn^('ic,i€)»  E<1uation  04). 


Page  1221  -  In  the  definition  of  a  , 
side  of  the  equation  should  be:  ! 


2im  N 


on 


l 

r=0, 1 


Equation  (36),  the  right-hand 


(C3) 


Part  II 


Page  1225  -  The  last  inner  bracket  of  Equation  (11)  should  involve 
the  derivation  of  the  radial  function  and  be  written  as  follows: 


(D* 


In 


(-ic,  i4 )  - 


R:  (-ic,io) 

in  R(3)(-ic»U) 

In 


Rj^(-ic.io) 


(C4) 


( I  \  *  MV 

Ron  (~ic»i°)  should  be  replaced  by  S^n'  (-ic.o)  as  noted 


Page  1226 

by  Flammer  in  the  errata  to  his  original  report  of  this  work. 


C.  Errata  in  Flammer's  Original  Report  [12] 


Page  12  -  Replace  sirs  by  coss  in  Equations  (49). 


APPENDIX  D 
COMPUTER  PROGRAM 


This  appendix  presents  the  Fortran  program  used  in  Section  VC 
of  this  work.  This  program  has  been  written  to  compare  Flammer's 
and  Meixner's  scattered  bistatic  far-field  in  the  normal  incidence 
perpendicular  polarization  case.  It  is  based  on  Hodge's  program[9] 
for  Meixner's  solution,  which  provides  the  structure  and  most  of  the 
subroutines  used.  The  reader  is  referred  to  reference  [  90  for  more 
complete  details  concerning  the  input  and  output,  and  for  structural 
questions.  First,  the  use  of  the  program  will  be  described.  Some 
remarks  will  then  be  made  on  the  programming,  and  a  listing  of  the 
program  will  be  provided. 

1 .  Use  of  the  Program 

The  program  uses  the  same  input  procedures  as  Hodge's  program 
except  that  the  two  following  data,  angle  of  incidence  and  polari¬ 
zation,  have  been  suppressed  since  only  the  normal  incidence  perpen¬ 
dicular  polarization  case  is  considered.  The  inputs  are  made  free- 
forrnat  on  a  teletype  as  follows: 

1.  KA  =  c,  electrical  circumference  of  the  disk. 

2.  THETA  SCATTERED  =  es  [degrees]. 

3.  PHI  SCATTERED  -  $s  [degrees]. 

As  in  Hodge's  program,  any  of  these  variables  can  be  incremented. 

The  corresponding  inputs  and  the  output  formats  are  identical  to  Hodge's. 
The  command  "ESC"  has  been  maintained.  It  interrupts  the  current 
calculation  and  requests  a  new  set  of  data.  No  negative  size  of  the 
disk  should  be  entered  and  a  disk  of  size  0  will  terminate  the  program. 

Unlike  the  general  case  treated  by  Hodge,  this  expresison  of 
the  far-scattered  field  in  the  normal  incidence  case  is  here  a  double 
infinite  summation  over  n  and  r,  the  index  m  taking  only  the  values 
0  or  1.  The  summations  are  truncated  in  the  same  way  as  in  Hodge’s 
work.  The  function  involved  in  eacb-term  of  the  summation  are  compared 
to  a  limit  fixedgin  the  program,  l(ru  in  this  case,  near  the  overflow 
level,  i.e.,  l(r°  for  the  computer  of  the  ElectroScience  Laboratory. 

To  prevent  any  overflow,  the  summation  is  truncated  when  one  of  the 
tested  functions  passes  this  limit.  This  procedure  allows  a  full 
use  of  the  capacity  of  the  computer  and  suppresses  the  need  for  a 
manual  search  for  the  optimum  truncation.  Flammer's  solution  involves 


the  spheroidal  functions  for  n-m  even  and  odd  and  their  derivatives. 

The  subroutines  computing  the  eigenvalue  of  the  spheroidal  scalar 
wave  equation  \  ( - ic)  and  the  expansion  coefficients  crn(-ic)  have 

been  extended  tonthe  n-m  odd  case  from  a  previous  work  of  Dr.  Hodge. 

The  subroutines  computing  the  functions  for  n  odd  have  been  added. 

2.  Remarks 

The  following  remarks  refer  to  the  program  listed  in  Section 
3  of  this  Appendix.  As  in  Hodge's  work  we  will  refer  to  a  line  of 
the  program  by  LN.  First,  the  different  parts  of  the  program  are  out¬ 
lined.  The  structure  is  similar  to  that  of  Hodge's  program. 

LN  10-53  Input 

LN  54-75  Computations  of  the  needed  functions 

LN  76-101  Computations  of  the  components  of  the  E-field: 

-  cross-section 

-  phase 

LN  102-104  Output 

LN  105-108  Incrementation  of  the  chosen  variable. 

The  needed  subroutines  are  listed  in  LN  113-564.  The  notations 
of  the  spheroidal  functions  for  n-m  even  as  defined  in  Hodge's  program 
[9]  have  been  preserved.  The  derivatives  of  the  functions  are  char¬ 
acterized  by  a  letter  "  P'  at  the  end  of  the  name  of  the  function  while 
a  suffix  OD  refers  to  the  function  for  n-m  odd.  The  following  notations 
must  be  introduced: 

A  and  B  are  the  two  summations  involved  in  the  computation 
of  and  b*. 

i 

ALPHA  and  BETA  in  the  main  program  are  equal  to  3  and  6  ,  re¬ 
spectively.  Those  names  come  from  Flammer's  definition  {7]. 

E0  and  El  are  the  summations  in  Equation  (112)  apart  from  the 
^-dependence  and  some  trigonometric  functions  of  0.  The  coef¬ 
ficients  necessary  for  the  computation  of  E0,  and  El,  are  com¬ 
puted  in  the  and  Rp^(-ic,io)  subroutines,  respectively, 

under  the  name°C0F.  n 


3.  Listinc 


In  this  section  is  given  the  listing  of  the  program  described 
above. 


1  I'PTjOrS  ,Vt« 

.•  r  F  At  FlH.I  t,r  f.TTE-Uf'i  BY  a  Cl»CU|  aK  FT  T  Al  L IC  n  I  RK 

f  i  t'iC I L)f *jf;t 

4  r  I  7  '  \,,j  c'.ll  A^  h ol  AUlSATlON 

ilMfi  uni  r<s.si-.SfS9 
t  II-- 1  'JUT  i>0  L7>'f  L(b)  ,VAR<T> 

7  CU'-PU  '<  f  4  iMPifUii  1  «  T  V  »  a  » t-  *  fOp  «  E  1  hr  ,  f.PH  I  ,  Ai.oHA  ,Bf  Ta 

(•  r0,'\'.)  t  t  It,  |  3  *i n  )  «l.  «  ’  'I a  .'<0  J  I'jFTb  l SO ) «F4  <  5n ) ,F4ptsO)  .COE  <  50  > 

9  1 .  bn  ( f  1J )  ,  S0(  (  Sf,  3  .P  <  50  )  ,  pH  (  bu  )  •  i-tT  A  (  50  )  iSt  TAOO  t  5n  ) 

IL  r>ATrt  L.'.litL/cHH  * THF.  R  PHT  5  / 

11  IA-(0..t.) 

IP  V.KT  rfc  ?*•  »'*<>) 

i pb  E  CO' •'.  'W  f  3  **///.*<  T  Tp£  »f  sr»  Tr>  BFSTAPT  PP.GORAw)../, 
l‘»  IMTYI’I  K/>  =  0  TC  STUP  PPOC.QaM)*) 

15  r  A|  L  F  .-r  (  V4P  ) 

lb  4  ror=i 

17  ICt_l  =  ll 

In  Inoexsl 

1°  l  K  I  T  f  ( >\*  »  «*  7  1 

,?0  ?1  v  i^rv  T<  ix«///,ix»  *1.  KA*«liiy«*s  •> 

21  l*LAO(B»->  vFP(l) 

PP  IMUAitin.LF.O.lFO  TO  5 

22  WHITE  <  <i  * «ic»  > 

ph  ?e  fui.iv.a  r  (u ,  tp.  theta  scattered  =  »> 

,'S  F!EA&to»->  VAM2> 

Pt-  lf.KTTKC'SP'3) 

27  p9  FOsmiATUX  »»••*.  PF I  SCATTfhEP  =  *) 

Pn  «L.',L'(5«“)  VAH(3» 

PC  HrtTTF <o,o0> 

in  ?o  ti>  ,//.!*, 'Which  variance  ts  to  »r  incremented?*) 

il  RErO(  :>«-»NVAR 

3P  IF  Ml\p/AR.Lt  .0  )  .OR*  CiVAP.liT  .(,)  )tiO  TO  ?3 

26  W*TTE<B,.U> 

-4  -.1  FUR'aaT  (l*t*TYPK  NUMBER  Op  CftcES:*) 

55  PEaC(c*->  ROC 

WKtTE»0,J2) 

37  ?«i  FOb'-iA  I  ( 1A  ,  •'ahaT  IS  the  ImcPEF’ENT?'  ) 

2?  RE  Ai  ,•  (  C  «  -  )  0  1 1'.'CRE 

’.V  Wi\VAK=  IVAh 

4;,  ILt  l=P*Hf"VAH-l 

**i  ?3  IUlP=TLI.l  +  l 

4/-  wk  r  TE  l  2  •*:**)  '-7  RFC  I  J  l  LI  )  ,  L/ibt  L  (  I  EL?) 

43  pH  rat'i'iA  r(i  7»//«px*2a- »lx« ’rKOSR  brcTlC'iv*  *?i  x*  *r  nORNS/i 

44  115y  ,  *  5tf,i'i ;•./  ( pi*A*  *>  )  *  ,!.ly  ,  *n  E  1  A  *  «  15X  *  *i  M  t  •  .  /  *  1 3X  * 

45  *PmI  *  *OXt  ’HAP*  .t>X.  *PHaSE*  **  X,  *,.1AG*  ,sX, 

hr  x  *  Phase •♦/ ) 

4 /  b  C=VAR ( 1 ) 

4b  «hi»,AX=45 

40  H|\l*  a X  —4 'j 

50  )  KI*VIAX  =  ';5 

M  Th-  S  =  v.  HPl*a.lH159/lno 

52  r I A=CbS( lut E» 

5  J  FhTSti/AM  (  j)  » :■*  .  T  HI  50/ lgo 

5H  (>U  10  H4=l,p 

55  F’sn-M-l 
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tib  lu-o 

5?  CAi  (.  O^Cibla  (C  .M»Hf'’'vAy  ) 

St  cau  u^rtf-i" ( c  A>. tr«N«.fiX,  tpkmax j 

V'C-  IF  IN) 

faO  CAl.t  bu,ijlF’i  f  ,'^rnAx  ) 

ht  Call  oL-f/.ut  (c a ,  j'i.mff'ax  .amKax,  i^r,waX  , 

ft,';  CALL  K!  L*f  (i  TA*M*]kK.iaX) 

l>.*  CA(  UI'Mm;,  { l;UK’A3(  *  1P9|>«AX) 

b«  CAlL  c  1  H.(m  (  .c o) 

bU  GO  ro  -7 

bf  CAI  L  ) 

A7  CAI  l  ULE«~N(l'tM,NI''i»l/.)0 

r,n  Cam  u  'hi r(c  ,r ,J,  H ,  (“mva*  ,LN«»AXt  tRRma*) 

by  Ca'  L  POLT.-.CETa*-"*  I^Ri-aX) 

7 u  CAI  l  OAAN(,  If  fjMAX.XPfipiAX) 

71  CAf  I  F  1 1  L|jJ  (  NF'N/FX  *  El  ) 

7i-  A 7  CONTINUE 

7  4  1  r,  CONTINUE 

7“  Ac»HAsiV(B-A  ) 

7  J  Pc.TA  =  -a/ ( Ff  -  A  ) 

7b  f  M;l  =  ALPIiA»I  0*SIMP.FXS> 

77  E'H^I=ALPHA4LC*rTA*C0S(PHrs>-[3E-TA*El*sIN/THEs>*r0s<PHT<;) 

7A  a  EhAi;TaCAt<S<ETMP) 

7^>  Lri/>oP=C  APS<  f  PhI  > 

Uu  SIi;thL  =  (  tiv  /f>T»2/C )  **<> 

fcl  Slr,pHl  =  (t.P'fc(  Pjitb/O***; 

a?  EKsHEaL ( tl PE  ) 

F l=rtIMAG(tTHf ) 

C‘4  IF  CL ff.KP.C.  >  GO  TO  ifc 

fib  AfcGsEX/Eh 

fib  FPf' AT  =  1  Aii/4  .  xul  59*  A  TAN  |  ApG  ) 

»7  iFMLK.Ll  .(  ,  )  .A'fl1.  If.  I.GT.u,  I  )  EPn/iT=FPHAT4.180 

»(■  IM  It  K  .  L  T  .  I, .  )  .A  Mi.  •  <Ll*l  T.o.  )  )EHHcTst  PbAT-lflO 

f>9  ifa  f>r/iTs:»o 

9»  IF  f  K'l  •!- T.n  .  >LPHATa-yo ' 

91  i  7  F-KrFfE/iL  I t.F'F'  t ) 

9P>  F  I  r:  A 1 1’1  A f}  ( t  F H  X  ) 

9.->  IF  <  LF  ,1'U.li,  J  GO  TO  ift 

9'»  AKpsEX/FK 

9b  EPpaF’s]  fi(//2. .  J  i*  t59*/‘  TAfj(  Ar&) 

9b  IFf  f  LK.Ll  ,(  .  > ,  AA'P.  IfTI  .GT.u.  )  > EpHAP=rPPAP4.l80 

97  jFf  (n:»Ll  •  A.  )  •  ANI/»  < <C J  .1  T.O.  )  jEpwtP=FpnAp»18n 

9o  GO  TO  19 

99  l  B  FPliAP=9o 

1  «><>  IF  f  cl  .LT.O.  )tPKAP  =  -90 

101  19  Tp  ( !JV A!.'  .  t  .  I  ’ )  ft’NVAKsl 

lo*  lkitei •A.i'M  vr.ntW'VrK)  .sJaxur  .sibf.nl  ,f  >V|A&r.ep«AT* 

1L-4  IF  M/WiP.FPHAP 

10A  P5  r-Of.  viat  I  jx.f  7.?,acly,Eir.*)  .2flA,EHJ.3«lv,p7.»)  1 

10b  3 7  IF  I  INUt  X.Ei».N0r>l'O  TO  4 

10F  J  NPF.Xsl  fil’i  y*l 

10  7  VAP  |  Ni«VA*  JsVAP  (NNVflO  I+VImcRE 

10F  GO  TO  « 

109  *>  C A|  l.  EX  1 1 

ill’  U  ^  CONTINUE 
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11 1 
11.' 

1 1  •>  r 
1  i  <•  r 
U‘>  r 
1  U  r 
117  r 
11'-. 

1 1-< 
lpn 
121 
12> 
12> 
124 
1.*  J  ? 
12.,  1 
127 
1  2.0 
12o 
l.i  r 
1JI 
l .I.-  r 

1J.)  r 
lJ4  r 

is  -j  r 

17,'..  r 

177 

lit 

IS-) 

lHo 

1M 

142 

14  7 
144 
14  j 
14b 
14  7 
14c! 

1  4  9 
1  j  • 

141 

142  r 
lJl  c 

1  J‘.  f 

161 
i-:.b 
IV.  7 
1  =  .-  4 

1.1') 
ll-r! 

1.1 
1  r»  •* 

1..7 

lu4 

It'S 


C-0  TO  4 
F  Ml 

ohi  aTk  .sf  hf  t-riPAL  A  i.jul.r"  f- umc  i  rr.c ,  s,  op  afagmEMT  n» 

OK.ipfi  ,  ■-ITM  M-)  fVFMl  or  Tu  CfiUFTi 

It-.  jAL  TO  t-’f  (0)) 

SU.  ^out  ji. c  '."f  n  (H,M.'ir,AX  > 
couple.::  F>,i4p,cOF 

rui  J'ry-  Flb(  1  r.  0  )  ,i,<  1  on  ,50  J  ,tir  rG(so)  »F4(5(;)  ,F4p(«cC>)  ,CnF<  =0) 

1  »  oil  (  5)U  ) 

SO ( 1 )=l 

T 1  jo  1 

nu  ,* 

SO  (  1  )  =  (  2*i"|-  -1  )*S»(  1  ) 
ru  .<  m.-j  , r.-j.-iAX 
I  •  =  !>  ♦  c  UI !  — A  ) 

V  U  (  l"'J  +  l  »  —  —  I  l.  +  M  +  l  )  *  SO  <  l|t!)  /  I  A'-|W  +  2  ) 

i  L'OC  [Iiuilf 

P  t.  T I  j  cl  <1 
Finn 

DtPIVATIVf.  (p  Dili  /,Tt  SF-NrhOlpAi.  aMGULaR 
FO-ir.TIOt  i  OP  AHbO..,EM  0;  OF-^f.ll  f.- ,  N , 

'■'Iti,  i-J-N  OuOS 

SO,. o oUT I>, e  Sf,MOtJ(i»ur'jf.,,viAX| 

f-O-  ,1,00  f  I  ,  I  Uitl  *(•(  1O0  ,'.0  j  ,C'rrG(so)  «F4(5r  )  ,F4p,rO)  ,COFi50  J 
1 ,Sr  ( so )  ,  GO)- ( bo  > 
c0 ■•.PLc.X  F  4  ,  F4p  ♦  cOF 
SU( (1)31 
FUrr-l  +  l 

(0  2  rfH  =  l,ftn 

SU(  ( 1  )  s  (  e  *|v  *,«1 )  *$0p  <  1 ) 

?  C  O’.  Tf  out: 

TO  7  MW=1 ,NM»#X 

sop  ( C;  j*l )  =-  I  rv  +  n*i.  >  »$CP  <  Mr  j )  /  t  M-1',  +  1 ) 

CO-  T  I MJF 

HCTUHI'1 

ENp 

Obi  aTE  S^hFROIDAL  EIGENVALUES  OF  flKGlif-'E.NT  C,  OROfR  M,n 

SO,  nfiOTir.E  Ot  EISM»C«Fi,A|NMaX) 

C Or.  .|0,’j  k'H.(Un) 

Ul«c;MSlON  IP(50) ,P(6G) ,A|  PhAfbO)  ,btTM50) 

coi,  rioi'F 

c«c=.r  >c 
Acr  =1 . nc-fb 
fJN'.-Miif.AXri 
i\  1  =>.||42  +  1 

P«1 )=1 
IH(1)=1 
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16  6 

IS-ft 

167 

1  C  Uf'/T  I  if’JK 

lob 

00  2  1)0=1,  i"N2 

16',. 

I  V-V»*  1‘JLi+.15>»1 

170 

lrtrM2  +  6*l(<^+lS 

17  J 

Ix  =  6'2  +  4*i,  ,  +2*ts*1 

17c 

Alp*  A  <  !  ia  )  =  (C?**1'**  (2*IV-i  ( IV-1  )-1  )  j/f  IX*  ( Ix-4  )  ) 

1 7  A 

i-  ( +  X  V-'l  1  M  )  V*'?  • 

17h 

? 

ntT  A  (  I‘jO+1  )  pr?/TX*f,(iHT  (lv*(lv  +  lr*lk»*(lu-l)/(tx*lx-'*.0)l 

17'_ 

f't-T  A  ( +  l  )  =0  . 

1;6 

BO- A'1$<  ALPf.Af  1  )  >+AFSIRFTfl<2)  1 

177 

00  6  Io'J=2,Vt!2 

175 

Atir-ABSft  LTA  (  I00)  )  *7  BS  j  ALphA  f  *04,  )  )  +AP?  (  Br  T  A  ( I  (3(3  +  1 >  ) 

17v 

PtT.il  1  -.."l )  -  ft  J  ft  «  Hi'J  )  *t‘CTA  1  1061 

If;  0 

IF  j  Ari.Or.bli)  BO  =  AU 

ltll 

A 

COMTlifUf; 

102 

AO=-ao 

1«  o 

13  o  T  =  B  0 

1C>4 

15 

C0MT1MUE 

Ifaj 

f’Ornol 

186 

00  2 0  Itfcsl  t.NMMAX 

187 

1=2* I jO-1+IS 

186 

N=1  +M-1 

186 

a=ao 

i9n 

B=PO 

191 

I tWRr •! 

192 

?1 

Ils=0 

19? 

Cgr(A+H)/2 

1?4 

IF'f  C0)1>0*22.50 

195 

50 

1 Kn  =  <0-A)/AiiS|C0) 

19c 

ItPRsiKKH+i 

197 

IF  (  If  K'i-ofi  )  '*C  ,'*1  *U1 

196 

Ul 

;vlIrTf  <  'l  ,**.?>  r ■ 

199 

02 

FOI.maTIjX,  »  ITfPATK.MS  EXrttt'FD  rCR  nOtCVALUr  •  ,  13* 

2  0  0 

C-0  TO  7ou 

201 

00 

lF(r9K“ACCl 24, 24,22 

2 IV 

?2 

B(P  )=AL  PI. A  (  1  )-C0 

20  A 

ro  5  1=8*M 

2C* 

n*  1 1  =  <  Af.no  a  <  i-i  )-co-bftA( 

205 

l  rtAf,  =  AOS<P<  I  )  1 

2  06 

IF  »PMAi..C,-r.)  .(If*35)li0  TO  7 

20  7 

5 

CO'fTINUF. 

20b 

J2 

i:OfTIiNi"F 

209 

no  6  1=2 «NI 

21C 

211 

6 

TFtPfl-n)‘'»9,14 

21c 

1* 

ItM  I  ) =*1 

2l3 

GO  TO  10 

214 

o 

iH(  n=i 

215 

10 

IF<lfJ(I1”lF’<  1-1  1)6,11,6 

216 

1* 

I i^sIIS+l 

21? 

6 

CONTIM'F 

2lt> 

IF  I  ns.-X«.t>lt.lb»l5 

219 

Asrn 

22  o 

GO  TO  21 
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221 

1  ^ 

f'=C0 

22; 

HI  TO  21 

2?  3 

?*» 

i’u=co 

2 24 

700 

tl:.(L)=-CO 

2^1 

20 

fOr>  T  I MJE 

22c 

JF(TS)  rtUftCOOiOUX 

227 

•  <i  n 

15  =  1 

22t- 

GO  TO  1 

22v 

nl  l 

CO'IT  IU'JK 

25< 

KltUR"1 

23; 

7 

M\tvAX=l-‘* 

23F 

< 1 1  =  PONm  AX  +  5 

23* 

f-0  TO  4 

2  34 

3  N|> 

23S, 

r 

*5l. 

r 

G3|  a  Tt  S^Hf'OinAL  t  l6Ef  FiiNCTTCim  fX^ftNSIOM  COrFP  1C  I  ENTS 

237 

r 

CF  Al<6'L)rf.NT  Ct  OHUtH  M.N.R, 

2  3:'. 

c 

23V 

FU'M.  OUT  TI'L  ChfPFMc.f  fM»1^AX,rN(<iAX*  IfipFAX  J 

24  0 

CO-.,»0,l  tlt-(loo)  *['<100,50) 

24  L 

r  I'xwiioN  ijF  (tom 

24  2 

C2-C*C 

2  43 

>1*1=  1+1 

244 

NN»’fl  ♦  ("IVF'Ax 

24b 

t  o  1  I\|N-1  ,Nf  MAX2 

246 

N  =  n+I'JN-1 

247 

J5=  (-1  )  **Ai|\J 

24  6 

IK(TS.I-E.O)  is=o 

249 

S  =  1S 

250 

4 

t'P(  TRK’AX+jlcO 

251 

f>(  jRio  AX  +  2)=1 ,«E-?0 

2  52 

f’  (  Mfj,  1  )  =0 

255 

L’lf  M*2),=l 

254 

00=<U'vi-l51/i  +  l 

255 

f'U  1.07  LU=l * IRRHAX 

2  5c 

<•  L  ”  1 

257 

IF  (1  L.(-r.OJ)  1  =IRR''1AX+jj.lL 

2  5r, 

Ift=>*L+l5 

2  5 '} 

i»\.-=v,+m 

2Co 

A«=<W+lR!-i  +  2)*|f!*lK''+l»*C?/|  |?*ihf+4»*<2  +  Ihm  +  «!)  J 

26  l 

f'K-  (2*X(>r'  »  (  J  4H  +  1 »  -2*i  *f>-l  )  +c?/<  |2*1Fw.*1>* 

262 

j 

.  * ?*  i„w+5 ) ) -rr*1*  <  i'o  +ti 

263 

CKsJH  +  1  1 4-1  )*C?S {  <  2* IRM-3) *  f  ?* IRN-i  J  ) 

264 

IF  ILU-  )j)n.t,10b,106 

265 

1  (.0 

rl ( i-., r  1 , t  +  3  )  s-f  rK»n<  A|*  .4  +3  )  +  ibR  +  EifijAif.  >  )*0f  r«.M,l  ♦  ?>  )/AR 

266 

Cl3A6=Al-:Fcr<!"M,l  +51  ) 

267 

IF  ( oft7*G-  tv_.  1  .  i.  .OF +30  )  jo  jo  3 

2al> 

60  TO  1(1 7 

2(,m 

106 

FF*  1 1  +1 )  =-  (  A,.  »pr>  <L  +  3  )  +  (PR+i  JOf  Nn  )  )  »DP ( l  +  ? )  )/Cp 

27, 

M'WsuiJSti.l  (  L  ♦  1  >  » 

271 

IF  oil, >1  .1 ,0f  ♦30) to  TO  3 

272 

1  C7 

COnT  I  6"  T 

275 

rir/ns (0(Ar..j.i+i» » 

274 

t’U=AL04lU<l’t_> 

275 

('LPjA’Tj  ( lip  (  6J  +  1  )  ) 
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27t  DUP  =  ft«-n&».f(l>V.P» 

277  ''L;4I>  jt'TL) 

t* 7 h,  nct'  r  Ai.i  5  ( 1-M.P  I 

279  rc=oi>"Uf 

2A0  IF  f  i)L  . r •  30  •  )  GO  To  s> 

2.  1  r  U:'=)  <  ,  du  f  1.  )  <  JJ+1 ) 

2d.-  f.L0  vi=«1*S(Ct;r.  ) 

*Cw  if- <  ACO'i.l’r  .1  ,cF  +  32)G0  TO  2 

2o9  TO  lit!  JS.Ju«  IPPMi>X 

2tV-  1  1  r  I  i..!(1.j*=?)=cr,M-4-f'c’Cj+2) 

2flS  F  =  1 

207  IS.c=IS>  +  N| 

anti  ik  r.  D  iyp.i9e»iy» 

2 r« o  1 9<i  f=3.o*m-i5> 

25,;  DU  1.1.0  1=1, iSM 

251  H’  F=F  *  t  IS!-,*  1 ) 


1  9(-  SU  i  =  0 

253  F'hV  =  IKRnAX  +  l 

291.  CO  no  1  =  1.  •‘■"IX 

29b  1K=P*X+Td 

29G  Syfi— So*l'*F  *  0  (  M9 ,  I* !•  > 

29/  TK(I-J>I)  nr-  ,157,133 

290  197  FlUf’sP 

299  113  F={ -F* ( 1K+2*M+1S-1) >/ ( JR»IS) 

300  ALPsPN^/SOFi 

301  CO  119  I  =  1,I'!9IX 

302  I'  (m.,4,  n=Ai_F*'.HNN»I  +  l) 

3 ;i 3  t  X  u  COi-'T  Itmlg 

3,19  1  CWMTIWJE 

4nb  PtrulVi 

30r>  3  lhP|V,AX=LL-l 

307  GO  TO  9 

30  f!  ?  1KP.«|AX=IHP|V'AX-1 

309  GO  TO  4 

31f)  b  NNvaX=«M|M-1  )/? 

311  HtTUKf) 

312  Fill! 

313  f 

319  r  MtfiATIWC  0  CUFFFICTEImT  SlJhRC'UTXMf 

31b  r  f--P.  iVF.N* 

3lij  C 

317  SUt.noUTIlvt  i.ifJtOWiC.f.MMMnX) 

316  (0«M0'-<  1.  lt><  lon>  »ul  j  OnifcO  1  ,  UNFG{5n> 

319  CO  u  i-iw=1  ,Nr.,v,flX 

32c  C2=C*C 

321  fi'l\  =  2*UM-l 

322  It-  (M,Gf  .11  GO  TO  ? 

323  OO  b  AiNsl,M9A,AX 

329  Ah=2*l\)'i-1 

3?b  $  tNri,(fM'-'»=U(FK,H 

32b  f-0  TO  3 

327  a  I’Asi.U 

32m  92=0.0 

329  Fi=flw«fJK» 

330  I'O  1  IPRsi.M 
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331  IK->*IRR-2-*l" -? 

34,'  A*-(2*  >ir  -*■:?)*  (  2*"i  +  lt>  +  1  >  *(  2/ <(  2*,, +2*  TrCPJS) 

344  i.*  (;.»•!♦>"*  ) ) 

334  t  r'  +  1  R  )  *  c  •  ♦r0*3  *  -K  1  -  C  ?*  C  r-  *rR  >  *1  +  1 ) 

i,i‘  i*cr>/<  ( **x  ,+<:*  i  r-i )  *  i  2*«-, +  2*  n-  ♦3  >  > 

3  3t>  ck-  ( in )  *  ( jR-i  j  *C2/  ( ( i:*r'*2*  i  r_4  )  *(2*i''+2*tk-1>) 

447  BOriC'? 

330 

43?  i  Px=(BR*[ia-c:;*n3»/AR 

44p  A=n  (  uk,  .  1 )  /si 

341  ON<_K(NfiJ=A 

34;  r  US-All-^  (L-HH’(NN»  > 

344  IP  Cr.MlK.l_T  .  1  .1  F-3b)r,0  TO  fe 

344  4  f  OuTI.iOfc: 

44t-  ^  Ct-rijRo 

340  (■  I" U,i  'VXr.wM-1 

34  7  At.  r URN 

344  t'No 

34  rJ  r 

3SC  r  Otic  ATt.  ■'OtOAC  rACIAL  FUNCTION  «C4)  OF  Af.rtUMENT  Cl 

4‘JI  r  CKmfH  K.1.5  l-  I  TO  fi--'  tvppj  I.P  TO  DKOFR  I\.sV:*2*MWF'A)r.2. 

3f>^  r  Al.Ro  I,UHM/*L14ATIUN  FUImCT tO!V*  m. 

3'ij  C  TL^o  COFUirifP'T  CCF  FOR  C0P  dIJ)  ATION  OF  f  ]. 

3f)4  r 

OD'J  SU'.'RkU’  Orr  AU<  C  If>«  “IP-MAX  •  IPK-IaX) 

300  C 0«-  y,um  f  luC  1(HI )  ♦!'*  1UC1  «LC  )  « Of-FG  ( F>0 )  ♦  F4  (  5m  * F^o  <  *0  )  *  C()F  (  SO ) 

i:-,7  1  ti>',  cu.O)  .SCl'CbO) 

3b<>  COiPLt.v  iX«P4»F**»i::«  COF  •  F4P 

35v  J x= { n .0 , 1.0  ) 

3  tul  PFaC  =  1 

36l  rAC  =  l 

3fs  F  Ar 1 

301  C-I\.,  =  l 

304  IHx.tU.UI  (-0  TO  2(| 

300  r  rty,V|  =  i;*i 

300  00  14  "IMsR.MAX" 

347  IKsMM-l 

364  k,Ko  =  |P.*lr-l)*|2*l'*)*(’hO 

3t>‘>  FF4C  =  l^4r1«"-i»i*P’frAC 

37c  FAc^s(  '  » I M  -  l  )  *  (  2  ♦  I  *;  )  *F  AC  ? 

371  PA(;sIC*fAc 

372  IP  (FFAO  ,l-f..l,pr*l  O00  TO  4 

373  IF  (  F  A  C  ,  0  *■  .1  .  0  F ♦ 3  U ) GO  TO  4 

374  2r,  |<S(0. 0*0.0) 

37b  10  )7  '-'Msj  ,NP;MAX 

37fc>  MV=>*fJ"-l 

377  Fi?*(N"'-l)+r 

470  J-U'  =  0 

37F*  ('Hrt.NO 

3o'>  FNf’H,4  =  n. 

3F.1  CO  1H  '  ;<  =  i,IrrmAX 

3o«  IK  =  2*(:r-1) 

3A3  r.UPPsM'tl’tijh  .NF  > 

i«4  £Uv=Soil.fSllp.P 

36b  I  i1/>C*=n*'?i t  01  T>K  ,NK  )  I 


3R0 

a  7 

3  Q  ‘t  p 
3hh 
39t 
391 

39;;  i  6 
334 

*  04 

3‘)d 

39* 

39/ 

396 

399 

4  01* 

HOI 

h  ,* 

H  J  3 
H  0  4 
Ht6 

Hot, 

HO/ 

Hr»; 

H(l9 

HU; 

HU 

H 12  1«* 

H13 

HJ4 

H 1 H  i  7 
Hit, 

HV7  1 
Hi  f. 

H 1 9  * 

HP.r 
H2l  5 
H2J* 

H._>.J  4 
H34 

H23 

H^t. 

H27  r 
H2a  r 
Hj*>  r 
430  r 
H3i  r 
H32  C 
HAS 
hah 
H3s 
H3G 
H  37 
H3b 

H39  2-) 

HHO 


IF  (!^i.'*.LT.1 .0F-SUH.O  TO  1 
H\lf;i»“K=2*r»K  ‘  <H(NK  *'**  )  )**p/<  ?*I,< 4*2 *M+  1  ) 
t  l\|OUi-’  =  F  f'ljt +  pf*  ORI"  P 
»*»S  (  ,,*(  ) 

1  K  { (V. ;:*.■•  K  i  .1  .  ft t  +  A  l  )  r-0  T0  * 

«'K=(Ih  +  .**l  +  A.)/|IW*J..)»(TH*24lM  +  2.,/(lH+?.)tr.ft 
CUMTJ  Jt ip 

Fl=(-1  >**,(■  K-1  t*J-***-»F7C*C**P*u(KK  ,1  )  /  (  (P*M+l  )*s>Lpi.') 

>*v*f  '-1  >»<<‘*'''-l>*FuC*r**(F-l>/(*MAC2)*s.lH159 

l*tFPA(.**2/SU  >*«>*'.  /L'.t*  r.frjNj 

F  h  . = A  U  S  (  P  P  ) 

IKfKRp.Pt  .  1.0F*6l')r,O  TO  4 
FF  Af=  <  >»*..  +1  )  »fL(7AC/  f  ) 

IF  |FFftrti)(  ,1.0F  +  1/)NAHAX=mN 
F  4  —  r(  1  “  I  X  *  K  p 
A Dz  Al1 6  (  FiJCihP  I 
“'UALtM'J  11  {  AC  ) 

AH  =  A<»S»kC) 

A2  =  aLuIj1i,(  A4) 

A3=Al+A2 

IF  (  ,» A .  i>T .  An  .  )  t- 0  To  3 

A3  =  AHf-(|!l) 

Al  =  ALU'*in(Af>) 

Aa-A96#Al>^AfiR(Ai» 

IF  (  .  r  .  .Mi  .  )  1,0  TO  3 

FH|;„iJ  =  t/(Fl'OR''**KH) 

(  0**t  If4*  'FT 

F £►*♦,>*  (-1  )#«hm*FH  «  ^'F#)*sO(Ai^')*SUl'#l 

COr  (IJ?J  )=-;»!)  *R1*FH  (  Mfj )  **S| j|x. 

CONTINUE 
FtTUlO'l 
FNnKNP=0 
00  TO  2 

FH||jr|)S(0,,0.» 

FU  TO  16 
GK  =  n 

r'U  TO  Id 
MIHPaXsi”' 

10  =  1 

KtTUR'X 

EWn 

I'Ll  iVATIVtS  OF  the  kAJ'IAi  FO^'C  T  TONS  pi  AND  PH 
OF  ORt'ER  N=n,N,  F  (.Dim  Oc  appUnent  0, 

ALSO  h!UK*Al  I/ATU-'M  FUr.'f  Tto.'i*  n. 

/'L<;o  tofFF  J<  iri'iT  C  f.)F  FOP  COMMUTATION  OF  FO, 

SU'-ROOT  |IMf_  C  HR  A  OP  <CM'!tA.ia,»HhAx.NNpiAX,|nRMAx> 

c  Or  NON  F  I0(  llit) )  »*’(  1  00 1 FO  )  ,L  Hrpf60)  »FH  (59)  .FHpf«iO)  *CoF  |  50) 
i*Sp|!-:0),suh(50) 

CO,PLtX  1X,PU,FH.RUF,F4P, a,C0F 
IXrtn. 0,1,0) 

F  F  A  C  s  1 
A=(0. 0.0.0) 
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